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ABSTRA£7r 


Hiis  thesis  considers  the  use  of  spreadsheet  techniques  as  a  foundation  for  the 
deveiopment  of  logistics  decision  support  systems.  An  inventory  model  is  presented  to 
show  the  flexibility  of  spreadsheet  techniques  and  demonstrate  the  use  of  various 
graphical  interface  techniques.  We  used  several  time  series  models  to  display  the 
graphical  capabilities  of  spreadsheet  programs  for  decision  making.  Finally  a  [)ro*v.ss 
control  model  is  prea^ted  and  its  implications  Hiscviss<Hh  Sample  spreadsheet  coding  is 
provided  for  ail  models  presented  with  a  primary  emphasis  on  graphic  output  to  aid  the 
logistics  manager  in  decision  making. 
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1.  INTRODUCTION 


A-  BACKGROUND 

Spreadsheet  programs  were  probably  the  most  significant  applications 
developed  in  the  early  days  of  the  personal  computer  revolution.  Originally 
designed  to  imitate  and  mechanize  the  accountant's  manual  spreadsheet 
calculations;  spreadsheets  were  readily  adopted  and  rapidly  gained  credibility 
in  the  area  of  Financial  Management.  By  applying  the  power  of  computer 
calculation,  in  a  familiar  and  time  tested  format,  the  financial  manager  could 
analyze  problems  that  had  previously  been  prohibitively  labor  intensive. 
Eliminating  the  drudgery  of  manual  calculation  allowed  the  Ifinancial  manager 
or  accoimtant  to  concentrate  of  the  problem  at  hand  and  has  resulted  in  a  vast 
array  of  spreadsheet  applications  that  support  the  financial  manager  in  his 
decision  making. 

The  continual  development  and  improvement  of  spreadsheet  packages 
have  further  contributed  to  their  importance  in  ftiancial  decision  making. 
Currently  available  packages  offer  a  broad  spectrum  of  special  functions  and 
graphical  features  that  allow  increasingly  complex  analyses  to  be  performed. 
Three-dimensional  spreadsheet  packages,  just  now  becoming  available  on  the 
market,  hold  even  more  promise  for  the  future.  Realizing  the  power  of 
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spreadsheet  analysis  and  a  certain  similarity  of  issues  between  finance  and 
logistics,  it  is  logical  to  c.  ,sider  whether  spreadsheets  might  also  be  useful  in 
the  area  of  logistics  decision  making. 

The  simplicity  of  modem  spreadsheet  packages,  the  availability  of  both 
those  packages  and  extremely  powerful  personal  computers,  and  the  increasing 
sophistication  of  logistics  managers  all  suggest  that  an  important  opportunity 
is  at  hand.  Current  spreadsheet  programs  contain  a  variety  of  specialized 
functions  for  the  automatic  handling  of  statistical  distributions,  present  value 
analysis,  table  analysis,  and  graphical  output.  These  specialized  functions 
offer  the  logistician  the  tools  he  needs  to  analyze  problems  and  develop 
specialized  applications  to  help  him  make  the  decisions  that  he  faces  every  day. 


B.  OBJBCnVB 

Several  problems  will  be  modeled  using  a  spreadsheet  program  and  the 
resulting  output  wiD  be  presented  in  graphical  format  as  a  demonstration  of 
the  spreadsheets  ability  to  help  the  manager  visualize  the  problem  and  arrive 
at  a  decision. 

It  is  recognized  that  spreadsheet  programming  may  not  be  the  most 
efficient  way  to  model  these  individual  programs.  A  specialized  program, 
produced  either  in-house  or  purchased  from  outside  vendors,  would  likely  be 
far  superior  to  a  simple  spreadsheet  program  developed  by  an  individual 
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manager.  However  the  prodigious  amount  of  time  required  and  the  high  cost 
associated  with  professional  software  development  often  makes  it  impossible 
to  produce  such  packages  in  a  timely  manner  or  for  a  limited  market.  The 
result  is  that  a  logistics  manager  often  makes  decisions  without  any  form  of 
decision  support  system,  relying  instead  on  intuition,  experience  and  luck. 

Spreadsheet  programming,  however,  offers  the  opportunity  for  any 
manager  to  write  programs  that  can  help  him  make  the  decisions  that  confront 
him  on  a  daily  basis.  By  studying  several  simple  problems,  the  value  of 
graphical  spreadsheet  output  can  be  evaluated  as  a  Logistics  Decision  Support 
Tool. 


C.  SCOPE,  UMITATrONS  AND  ASSUMPTIONS 

The  primary  focus  of  this  thesis  is  on  the  potential  use  of  spreadsheet 
graphical  output  as  a  basis  for  logistics  decision  support  systems.  The  value 
of  both  graphic  and  tabular  output  are  considered,  and  various  examples  are 
provided  to  demonstrate  the  ability  of  the  spreadsheet  to  help  a  manager 
better  understand  a  problem  and  make  better  decisions. 

Microsoft  Excel  Version  4.0  was  selected  for  this  stu(ty  because  it  is  one 
of  the  newer  additions  to  the  list  of  avedlable  spreadsheets  on  the  market  and 
offers  a  variety  of  functions,  a  comprehensive  macro  language  and  extensive 
graphical  capabilities.  In  addition.  Excel  is  reasonably  easy  to  learn  and  offers 
a  variety  of  programs  to  translate  its  code  into  other  spreadsheet  formats. 
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Coding  was  not  considered  particularly  important  to  this  thesis  and  may 
not  reflect  optimal  or  even  good  coding.  The  primary  ijitent  was  to  develop 
code  that  worked  sufficiently  well  to  provide  graphical  output  for  consideration. 
Portions  of  the  code  included  in  the  appendices  are  provided  so  that  others 
more  interested  in  code  generation  can  have  a  starting  point  from  which  to 
consider  the  problems  presented.  The  problems  analyzed  in  the  body  of  this 
thesis  were  selected  to  illustrate  the  use  of  spreadsheet  graphical  output  and 
were  not  necessarily  selected  for  their  individual  value  in  logistics  decision 
making.  Probiems  were  primarily  drawn  fi:t)m  existing  literature  and  are  well 
documented.  The  data  used  were  generated  by  spreadsheet  models  using 
statistical  techniques  and  may  not  be  applicable  to  real  world  scenarios 
without  modification. 

D.  METHODOLOGY 

This  research  relied  heavily  on  existing  models  and  theories  relating 
primarily  to  confidence  interval  procedures,  process  control  procedures  and 
acceptance  sampling  procedures.  The  models  developed  were  adapted  firom 
those  presented  in  several  papers  and  implemented  as  spreadsheet  programs. 

Though  not  specifically  addressed,  the  power  of  current  spreadsheet 
programs  is  quite  apparent.  Much  of  the  work  on  which  this  paper  is  based 
was  originally  conducted  using  extensive  mainframe  computer  assets  and 
considerable  time  in  developing  suitable  output  formats.  Much  of  this  can  now 
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be  replicated,  with  minimal  time  and  effort,  on  a  personal  computer  using  a 
relatively  simple  spreadsheet  program.  The  graphical  capabilities  of 
spreadsheet  programs  allow  the  decision  maker  to  quickly  and  easily  display 
the  results  of  his  work  so  that  he  can  better  \dsualize  the  system  he  is 
studying. 

B.  ORGANIZATION  OF  THE  STUDY 

Chapter  II  discusses  the  background  of  the  problem.  In  Chapter  III  an 
inventory  simulation  model  was  developed  using  spreadsheet  programming. 
In  Chapter  IV  several  time  series  models  were  studied.  Chapter  V  discusses 
the  application  of  spreadsheets  for  statistical  process  control  with  correlated 
observations.  Chapter  VI  provides  concluding  remarks. 
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11.  BACKGROUND  OP  THE  PROBLEM 


Rapid  developments  in  the  field  of  logistics  have  lead  to  an  increasingly 
complex  environment  for  the  logistics  decision  maker.  Once  merely  a  loose 
collection  of  procedures  used  to  support  deployed  military  operations;  logistics 
is  now  developing  into  n  science  with  far  reaching  implications,  and 
applications  far  beyond  the  realm  of  military  support.  While  logisticians  were 
once  concerned  primarily  with  moving  material  they  are  now  faced  with 
deciding  what  material  to  move,  how  to  move  it,  the  likelihood  it  will  be 
required  if  it  is  moved  and  a  host  of  other  decisions. 

Because  of  this  increase  in  scope  the  logistician  is  now  faced  with  a 
perplexing  array  of  decisions,  many  of  which  are  no  longer  intuitive  or 
experience  based.  Some  form  of  decision  support  system  must  be  adopted  to 
help  in  making  timely,  efficient  and  correct  decisions.  No  longer  are  pencil  and 
paper  solutions  sufficient  to  solve  the  problems  presented.  The  time  required 
for  manual  calculations  could  alone  render  the  decision  worthless,  and  the 
potential  for  error  in  such  manual  calculation  could  be  even  worse. 

Mainframe  computers  and  calculators  offered  the  first  solutions  to  the 
problem  but  neither  was  well  suited  to  the  task.  Mainframes  were  too  big,  too 
inflexible,  and  too  complicated  for  use  in  other  than  system  wide  decision 
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making.  Calculators  were  portable  and  easy  to  use  but  not  sufficiently 
automated  or  powerful  when  dealing  with  complex  logistics  problems  and 
concepts  that  now  confront  individual  managers.  Personal  computers  offered 
the  first  realistic  hope  of  providing  logisticians  with  a  truly  useful  and 
universal  platform  for  the  development  of  decision  support  systems. 

Development  of  personal  computers  was  insufficient  by  itself  to  solve  the 
problem  however.  The  initie?  personal  computers  functioned  very  like  the 
mainframes  on  which  they  were  modeled.  They  were  slow,  not  readily 
available  and  generally  provided  cryptic  output  at  best.  It  was  the 
development  of  the  spreadsheet  program  that  ultimately  offered  true  hope. 

Spreadsheets  were  initially  modeled  after  manual  documents  used  by 
accountants.  Accountants  had  used  the  same  form  for  hundreds  of  years,  a 
ruled  sheet  of  paper  with  columns  and  lines  that  intersected  to  form  little 
rectangular  boxes  into  which  they  wrote  the  information  pertaining  to  the 
problem.  The  numbers  in  the  boxes  could  then  be  added  up,  either  by  column 
or  by  row  and  final  answers  could  be  figured  out.  By  automating  that  simple 
process  on  the  computer  untold  hours  of  calculation  and  copying  were 
eliminated.  Problems  that  were  previously  not  considered  worth  the  cost  of 
obtaining  a  solution  could  now  be  programmed  and  easily  handled. 

Financial  managers  immediately  adapted  to  the  new  technology.  The 
format  had  been  designed  for  them,  was  familiar  and  eliminated  the  drudgery 
of  calculation.  Logisticians  were  then  just  starting  to  develop  the  complicated 
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ideas  to  improve  the  performance  of  logistics  systems.  As  they  worked  to 
change  logistics  from  an  art  form  into  a  complex  science,  they  rapidly 
developed  a  need  for  automated  systems.  But  the  concepts  hadn't  been 
developed. 

Gradually  over  the  last  several  years  logisticians  developed  the  necessary 
scientific  structure.  People  were  trained  in  probability,  math,  calculus,  finance 
and  a  myriad  of  other  disciplines.  The  range  of  problems  claimed  by 
logisticians  as  within  their  domain  increased  to  the  point  where  complex 
specialization  was  essential  and  some  form  of  computer  assistance  was 
required.  Unforttmately  the  rest  of  the  society  didn't  realize  what  was  going 
on. 

Even  the  author,  with  twenty-one  years  of  experience  in  a  variety  of 
military  supply  and  logistics  positions,  didn't  know,  until  quite  recently,  that 
logistics  is  now  concerned  with  probabilities,  failure  rates,  availabilities,  and 
the  many  other  esoteric  ideas  logisticians  have  chosen  to  adopt.  Ordinary 
people  think  that  logistics  is  still  the  art  of  moving  material  to  support 
deployed  military  forces.  Consequently  no  universal  logistics  software  has 
been  forthcoming  to  solve  the  problems  that  confront  the  logistician.  Few 
people  outside  the  specialty  xinderstand  that  such  software  is  required,  after 
all  just  moving  material  doesn't  require  complex  computer  techniques. 

There  have  been  many  programs  developed  to  solve  specific  logistics 
problems,  but  they  typically  require  the  user  to  manipulate  the  data  into  a 
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form  acceptable  to  the  particular  program  and  only  provide  some  limited 
output  that  may  not  be  of  much  real  use  to  the  decision  maker.  Custom 
programs  can,  and  have  been,  written  for  various  logistics  problems  but  they 
have  traditionally  been  long  laborious  projects,  written  in  some  highly  complex 
computer  language  by  professional  programmers  who  have  little  understanding 
of  the  purpose  of  the  programs  they  write.  Little  has  truly  been  done  to 
provide  the  logistician  with  ready  access  to  a  broad  range  of  programs  essential 
to  his  work. 

Possibly  this  problem  has  finally  been  solved  by  the  new  generation  of 
spreadsheet  programs  currently  available.  No  longer  restricted  to  one  screen 
of  columns  that  can  be  acted  on  by  a  limited  number  of  functions;  the  modem 
spreadsheet  provides  a  broad  array  of  functions,  the  capability  of  having 
multiple  interrelated  spreadsheets,  a  simple  macro  language,  and  extensive 
graphics  capabilities.  The  modem  spreadsheet  may  be  the  software  that  offers 
the  logistician  the  power  exr^>it  the  personal  computer  just  as  financial 
managers  have  been  doing  for  years. 

The  logistician  now  has  the  tools  to  develop  complex  logistics  decision 
support  systems  on  his  own.  By  using  simple  spreadsheet  commands  the 
logistician  can  now  design  systems  to  support  him  in  his  decisions  as  the  need 
for  them  rises.  He  can  personally  take  control  of  the  personal  computer, 
without  the  interference  of  computer  programmers,  and  develop  programs  that 
actuaUy  help  him  make  decisions. 
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Perhaps  the  most  innovative  feature  of  this  new  generation  of  spreadsheet 
programs  is  the  advances  that  have  been  made  in  graphics.  Once  a  tedious 
process  to  obtain  a  barely  legible  graph;  one  can  now  by  using  a  few 
keystrokes,  or  a  mouse  interface,  create  clear,  legible  and  informative  graphic 
representations  vdth  ease.  These  can  be  particularly  important  in  analyzing 
complex  logistics  relationships,  for  decision  making  or  for  educational  purposes. 
The  remainder  of  this  paper  will  concentrate  on  the  potential  of  this  graphical 
output  as  a  foundation  for  Logistics  Decision  Support  Systems. 
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ni.  AN  INVENTORY  MODEL 


The  idea  of  using  spreadsheet  programming  for  financial  analysis  is  not 
new.  Procedures  have  been  identified  and  documented  as  guides  for  managers 
to  employ  in  developing  spreadsheet  programs  to  aid  them  in  financial 
analysis  and  decision  making.  In  a  paper  on  spreadsheet  risk  analysis  (Seila 
and  Banks,  1990),  a  detailed  example  is  offered  to  demonstrate  the  use  of 
simulation  techniques  to  develop  a  risk  analysis  model  from  an  existing 
spreadsheet  through  the  use  of  spreadsheet  macros. 

Simulation  is  a  useful  method  for  studying  real  world  systems  that  are 
either  too  complex  or  too  reliant  on  random  probability  for  analytical  analysis. 
Simulations,  employing  models  based  on  the  important  variables  involved, 
provide  considerable  insight  into  the  operation  of  such  systems  and  can  serve 
as  a  foundation  for  decision  making  when  more  precise  methods  are  not 
available.  While  simiilation  can  be  performed  manually,  the  complexity  of 
problems  that  may  be  studied  using  manual  techniques  is  extremely  limited. 
Spreadsheet  programs  and  macro  techniques  offer  the  manager  a  powerful  tool 
for  analyzing  such  problems. 

Spreadsheet  packages  offer  sophisticated  mathematical,  statistical  and 
financial  functions  which  the  manager  can  use  to  develop  models  of  more 
complex  systems.  Graphical  capabilities  can  provide  important  visual  aids  for 


the  decision  maker  as  he  seeks  to  identify  trends  or  patterns  inherent  in  such 
systems.  Database  functions  can  provide  the  manager  with  ways  to  catalog 
and  store  his  results,  and  macro  programming  can  automate  much  of  the 
process. 

Macros,  which  are  short  programs  written  to  accomplish  specific  tasks, 
can  be  used  to  automate  complex  procedures  and  are  ideally  suited  for  the 
execution  of  the  repetitive  calculations  typical  of  simulation.  Using  macro 
programming  a  manager  can  quickly  and  easily  analyze  a  number  of  competing 
scenarios  and  improve  his  understanding  of  a  problem.  The  spreadsheet  can 
thus  become  a  tool  to  help  the  manager  make  better  decisions  by  increasing 
knowledge  of  how  various  systems  are  likely  to  perform. 

An  inventory  model  was  constructed  to  explore  spreadsheet  simulation  of 

logistics  problems  using  Microsoft  Excel  Version  4.0.  This  problem  is  typical 

of  those  presented  in  logistics  classes  and  is  presented  as  follows; 

(Inventory  System  with  Stockouts  and  Backorders):  A  Navy  Supply 
Center  (NSC)  is  planning  to  install  a  system  to  control  the  inventory  of 
a  particular  nonrepairable  (consumable)  item  caUed  component  XYZ.  The 
time  between  demands  for  XYZ  is  uniformly  distributed  between  0.2 
weeks  and  0.8  weeks.  The  requisition  amount  of  each  customer  is 
uniformly  distributed  between  2  and  6.  In  the  case  where  custonlers 
demand  XYZ  when  it  is  not  in  stock  (a  Stockout),  80  percent  must 
arrange  a  special  order  at  100  dollars  per  order  whenever  the  demand 
occurs,  while  the  other  20  percent  are  not  urgent,  thus  will  backorder  and 
wait  for  the  next  shipment  arrival.  NSC  employs  a  periodic  review- 
reorder  point  inventory  system  where  the  inventory  status  is  reviewed 
every  four  weeks  to  decide  if  an  order  should  be  placed.  NSC's  policy  is 
to  order  up  to  the  stock  control  level  of  72  XYZ's  whenever  the  inventory 
position,  consisting  of  items  in  stock  plus  items  on  order  minris  the  items 
on  backorder,  is  found  to  be  less  than  or  equal  to  the  reorder  point  of  18. 
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The  procurement  lead  time  (the  time  from  placement  of  an  order  to  its 
receipt)  is  constant  and  requires  three  weeks.  The  initial  inventory 
values  are  72  units  on  hand,  0  units  on  order  and  0  units  on  outstanding 
backorder.  (Modifred  from  Pritsker,  1984,  pp.l94-196) 

Although  a  somewhat  simplified  version  of  reality,  this  problem  is 
illustrative  of  the  complexity  involved  in  typical  logistics  problems  for  which 
decision  support  systems  are  required.  There  is  obviously  no  single  answer 
that  will  provide  a  universal  solution-  Consequently  a  suitable  decision 
support  system  would  serve  primarily  to  provide  a  clear  understanding  of  the 
processes  involved  rather  than  a  single  answer.  Such  a  system  could  help  a 
manager  make  better  decisions  than  might  otherwise  result. 

A  simple  reading  of  the  problem  shows  that  there  are  a  number  of 
variables,  operating  together,  which  affect  the  inventory  at  any  point  in  time. 
By  analyzing  these  variables  and  the  likely  relationships  between  them  a 
model  of  the  system  that  closely  resembles  inventory  behavior  in  a  real  world 
situation  can  be  developed.  A  manager  can  then  use  a  spreadsheet  program 
to  code  these  relationships,  automate  the  calculations  involved,  and  design 
graphical  interfaces  to  help  him  understand  and  make  inventory  decisions. 
Code  developed  for  the  inventory  model  can  be  foimd  in  Appendix  A. 

The  primary  benefits  of  such  an  approach  are  the  virtual  elimination  of 
calculation  errors,  a  tremendous  reduction  in  the  time  required  to  perform  the 
calculations,  and  a  considerable  flexibility  in  displaying  results.  Additional 
benefits  include  easy  modification  of  the  model  as  new  variables  are  identified. 
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adaptability  of  the  model  to  similar  type  problems,  and  through  the  use  of 
macros  the  ability  to  provide  a  user-firiendly  system  for  less  knowledgeable 
personnel  to  use. 


One  of  the  features  offered  in  Excel  4.0  is  the  ability  to  generate  a 
graphical  interface  through  which  users  can  modify  specific  program  variables 
and  alter  program  operation  without  recoding  the  underlying  spreadsheets. 
This  feature,  called  the  dialog  box,  can  be  coded  into  a  macro  whenever 
exploration  of  various  options  might  be  of  interest  to  the  user.  This  greatly 
enhances  the  users  options  and  allows  him  to  concentrate  on  the  results  of  the 
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Figure  1:  Input  Dialog  Box,  Inventory  Model 
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model  rather  than  the  code.  The  dialog  box  developed  for  operator  input  in  the 
inventory  model  is  shown  in  Figure  1. 

As  seen  in  Figure  1.  the  user  has  an  option  to  change  distributions  as  well 
as  parameters  of  the  requisition  order  frequency,  the  requisition  order 
quantity,  stockout  information  and  procurement  lead  time.  The  Input  Sheet 
provides  two  types  of  inventory  system  for  analysis  both  the  R- System,  which 
bases  reorder  frequency  on  time  since  the  last  reorder,  and  a  Q-System  which 
bases  reordering  on  inventory  level  reaching  a  predetermined  level.  There  is 
also  an  option  to  change  the  initial  inventory  values,  time  period  to  record  and 
various  other  values.  Other  options  could  be  provided,  at  the  discretion  of  the 
individual  who  codes  the  spreadsheet. 

Using  a  dialog  box  as  an  interface  for  user  input  allows  the  spreadsheet 
programmer  an  ability  to  both  request  information  firom  the  user  and  to  verify 
the  format  of  the  users  response.  The  input  screen  used  in  the  inventory 
model  provides  twenty  three  options  which  can  be  varied  by  a  user  interested 
in  exploring  different  scenarios.  Various  options  include  modifying  the 
distribution  of  random  variables  generated  by  the  simulation,  setting  initial 
inventory  values,  and  specifying  the  type  and  nature  of  the  inventory 
replenishment  system  being  considered.  There  is  also  a  provision  for  obtaining 
screens  to  guide  the  user  in  entering  specific  values  and  a  cancel  option  to  stop 
further  program  execution.  A  similar  interface  was  designed  for  information 
output  (see  Figure  2). 
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Although  not  as  elaborate  as  the  Input  Dialog  Box  the  Output  Box 
provides  an  equally  important  interface.  The  spreadsheet  programmer  can 
offer  a  list  of  optional  outputs  which  the  program  (!an  provide.  The  user  can 
then,  by  selecting  a  button  on  the  screen,  obtain  any  of  the  available  »)utputs 
desired.  Using  dialog  boxes  in  this  way  allows  spreadsheet  programs  to  be 
developed  and  used  by  managers  with  varying  degrees  of  expertise  in 
spreadsheet  programming  and  logistics  problem  solving. 


Figure  2:  Output  Dialog  Box,  Inventory  Model 

The  coding  for  each  of  these  dialog  boxes  is  included  in  Appendix  A. 
While  it  may  appear  somewhat  complicated  it  is  easily  generated  using  a 
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separate  program  provided  with  Excel  4.0,  the  Dialog  Editor.  The  Dialog 
Editor  is  itself  a  graphical  interface  that  allows  the  spreadsheet  programmer 
to  design  a  dialog  box  on  screen,  using  drag  and  drop  procedures.  When  the 
screen  is  completed  it  is  saved  and  the  coding  required  is  generated  by  the 
Dialog  Editor.  This  code  can  then  be  called  by  a  macro  program  when 
necessary  to  activate  the  user  interface. 

Examples  of  the  graphic  output  are  provided  in  Figures  3  through  6. 
These  examples  show  graphs  that  might  be  generated  by  a  decision  maker  who 
was  interested  in  the  effect  on  inventory  levels,  backorders  and  periods  of 
inventory  stockout.  By  generating  graphs  using  two  different  assumptions 
concerning  the  level  to  which  inventory  must  drop  before  a  reorder  is  placed 
the  decision  maker  can  compare  various  scenarios  and  likely  improve  his 
decision  making.  Figures  3  and  4  graph  the  Inventory  Position  (On  Hand  +  On 
Order  -  Backorder)  and  Figures  6  and  6  graph  the  Net  Inventory  (On  Hand  - 
Backorder)  simulated  by  the  model  for  two  different  reorder  points,  of  18  and 
10  respectively.  All  other  variables  are  held  constant. 

We  expect  a  higher  inventory  level  when  the  reorder  point  is  18  than 
when  it  is  10,  we  also  expect  less  stockouts  with  the  higher  reorder  point. 
Comparing  Figures  5  and  6,  we  can  see  that  the  length  of  stockouts  and  the 
total  number  of  items  unavailable  when  requested  is  greater  with  the  lower 
reorder  point  (Figure  6).  Also  we  can  see  that  the  maximum  amount  of 
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inventory  on  hand  at  the  lower  reorder  point  (Figure  6)  is  typically  much 
lower. 

By  comparing  these  4  graphs  a  manager  can  observe  a  number  of 
important  things  about  how  his  inventory  will  react  to  various  reorder  points. 
Together  these  graphs  can  help  him  visualize  and  understand  how  the  reorder 
point  affects  the  number,  size  and  frequency  of  stockouts,  the  size  of  orders 
average  inventory  and  a  variety  of  other  detail.?..  These  graphs  don't  provide 
a  tlnai  solution,  but  they  can  be  instrumental  in  providing  the  manager  with 
the  insight  he  requires  to  make  a  decision. 

This  type  of  graphical  support  can  provide  the  decision  maker  with  an 
idea  of  what  is  going  on  in  the  system  he  is  simulating.  He  can  see  the  high 
and  low  inventory  positions  graphically  displayed  based  on  the  parameters  he 
has  specified  in  his  model  over  a  series  of  trials.  Further,  he  has  an  option  to 
change  the  parameters  and  to  nm  the  simulation  again  to  analyze  the  effects 
of  various  changes  on  the  high  or  low  inventory  positions  he  is  interested  in. 
This  information  can  help  him  iinderstand  the  impact  of  each  variable  on  his 
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manager 


is  better  able  to  make  decisions  related  to  this  type  of  system. 

The  decision  maker  Jilso  has  ready  access  a  the  same  data  in  tabular  form 
as  shown  in  Figure  7. 

Reorder  R=  18  Reorder  Pi=  10 


Average  Customer  Requisition  Size  =  3.6566  3.3667 

Number  of  Reorders  ■  B  B 

Average  Requisition  Size  =  55  5  56  333 

Number  of  Stockouts  =  33  71 

Average  Net  Inventory  -  13.098  1  75 


Figure  7;  jutpuc  raoie 

By  viewing  the  results  in  tabular  form  the  manager  can  consider  the 
specific  values  involved  in  greater  detail.  The  simulated  average  requisition 
size,  number  of  stockouts,  number  of  reorders,  etc.  can  all  be  easily  determined 
by  reviewing  the  tabular  output.  In  addition  this  format  output  format  may 
be  more  compatible  with  the  decision  making  style  of  some  managers.  From 
Figure  7,  it  is  clear  that  the  policy  of  setting  the  reorder  point  to  18.  as 
compared  to  a  reorder  point  of  10,  yields  a  smaller  number  of  stockouts  at  the 
expense  of  maintaining  a  higher  inventory  level. 

Spreadsheets  offer  the  ability  to  manipulate  data  into  forms  that  are  most 
useful  to  an  individual  decision  maker.  In  any  case  the  value  of  spreadsheet 
analysis  of  logistics  problems  is  significant.  A  complex  problem  can  be 
simulated  or  actual  data  finm  historical  observation  can  be  input  and  the 
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simulation  portions  of  the  model  bypassed  The  resulting  data  can  then  Ije 
analyzed  or  displayed  rapidly  by  the  manager  using  spreadsheet  commands  fo 
group,  organize  or  present  the  data  in  ways  that  can  provide  meaningful 
insight  and  help  improve  decisions  . 
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IV.  SEVERAL  TIME  SERIES  OP  INTEREST 


A.  INTRODUCTION 

We  will  now  consider  applying  spreadsheet  techniques  to  a  less 
straightforward  logistics  concept.  Statistical  analysis  often  demands  the 
assumption  of  independent,  identically  distributed  (i.i.d.)  random  variables. 
Often  this  is  not  the  case.  Logisticians  frequently  deal  with  such  variables  as 
units  produced  on  only  one  machine  or  c'omponent  failures  occurring  in  one 
particular  type  of  aircraft.  Are  these  events  independent?  In  fact  these  events 
may  be  correlated-  Either  the  machine  wears  at  a  constant  rate  or  the  aircraft 
characteristics  impose  a  particular  stress  on  a  given  part.  Any  of  a  myriad  of 
details  can  influence  the  random  nature  of  the  obsen’ed  part.  If  the  systematic 
error  involved  is  significant  the  rule  of  independence  is  violated  and  the  system 
under  consideration  may  not  be  amenable  to  a  simple  statistical  analysis. 

This  problem  has  been  frequently  discussed  in  the  literature,  but  except 
in  an  academic  environment  such  readings  are  seldom  high  on  a  practitioners' 
reading  list.  It  is  much  simpler  to  assume  that  all  variables  are  both 
independent  and  identically  distributed  and  more  restrictive  cases  are  seldom 
considered  at  the  practiced  level.  However,  by  applying  spreadsheet  technology 
to  this  notion  it  is  easy  to  see  problems  associated  with  this  simplification. 
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Using  ideas  developed  in  Kang  and  Schmeiser  (1990).  a  graphical  analysis  of 
several  correlated  processes  will  be  considered. 

We  will  use  three  different  time  .series  models  to  understand  the 
difficulties  in  analyzing  correlated  data:  AR(1)  model,  EAR/ 1 )  model  and  M/M/1 
queuing  model.  By  using  simple  spreadsheet  techniques  each  of  these 
situations  can  be  easily  modeled  and  the  results  graphically  displayed,  results 
which  will  immediately  convince  even  the  most  skeptical  observer  of  the 
potentially  undesirable  results  obtained  if  an  erroneous  assumption  of 
independence  is  employed. 

B.  AR(1)  MODEL 

The  AJR(1)  process,  can  be  defined  by  the  equation; 

Xt  =  ^  (1) 

where  -  N(0,l-(t)^.  Without  loss  of  generality,  we  set  g  =  0.  Note  that  AR(l) 
process  with  4)  =  0  becomes  normal  i.i.d.  process.  As  (j)  approaches  to  1  fi:t)m  0 
higher  correlation  will  occur.  The  steady  state  AR(1)  process  can  be  generated 
by  setting  Xj  ~  N(0,1)  and  recursively  generating  Xg,  X^, . 

Using  a  spreadsheet,  graphic  representations  of  such  a  system  can  be 
generated  to  provide  a  visual  presentation  of  the  typical  distribution  of  such 
a  process.  Figure  8  represents  a  series  of  100  sequential  values  generated 
using  Equation  (1),  with  (1)=0.9. 
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Figure  8:  100  Observations  from  the  AR(1)  Model  with  <j>=0  . 


A  100(1-4))%  Confidence  Interval  can  be  constructed,  using  the  equation 


y/n 


(2) 


where  X  is  the  sample  mean,  defined  as 


;  and  S  is  the 


sample  standard  deviation,  defined  as  5^ = 


t 


(Xj  -  JOV(^-i) ;  t,.^  „.iisthe 


(l-(V2)th  quantile  of  the  Student's-t  distribution  with  n-1  degrees  of  freedom; 
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and  n  is  the  number  of  observations  in  the  sample.  The  quantity  -S/ /n  is 


referred  to  as  the  standard  error. 

To  further  explore  the  properties  of  an  AR(1)  model  a  number  of  points 
are  generated  and  graphed  following  a  procedure  described  in  Kang  and 
Schmeiser  (1990).  A  spreadsheet  program  was  created  to  generate  4,096 
sequential  values.  These  values  are  then  grouped  into  successively  larger 
batches  with  the  resulting  mean  value  of  each  batch  being  considered  as 
representative  of  that  batch.  The  first  16  mean  values,  for  each  batch  size,  are 
then  used  to  represent  16  samples  of  that  given  batch  size  drawn  successively 
firom  an  ongoing  process. 

This  procedure  is  replicated  for  100  times.  The  confidence  interval  can 
be  calculated  using  Equation  (2).  and  the  coverage  for  a  given  nominal 
coverage  1-a  (we  use  a=0.1  as  an  example)  is  estimated  by  the  proportion  of 
the  points  that  lie  above  the  lines  (Kang  and  Schmeiser,  1990), 


where  p=0  in  this  case.  In  Figure  9,  100  replications  of  16  batches  each  of 
which  has  a  batch  size  1  (i.e.  m=l)  were  simulated  for  the  AR(1)  process  with 


(t)=0.9.  The  coverage  is  estimated  by  the  proportion  of  the  points  that  fall  above 
the  V-shape  lines  calculated  from  Equation  (3)  that  is  shown  in  the  figure.  The 
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graph  shows  that  the  values  are  widely  dispersed  and  that  only  21%  of  the 
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while  the  nominal 


coverage  is  90%.  As  mentioned  before,  this  pofir  coverage  is  due  to  the 
correlation  between  observations. 


Figure  9:  .AR(1)  (jnO.S  Batched  in  size  m=l,  with  16  batches 

When  the  process  contains  positive  correlation  as  can  be  observed  in  this 
case  of  AR(1)  model  with  4»  =  0.9,  S^,  defined  in  Equation  (2),  is  not  an 
xinbiased  estimator  of  o^,  the  process  variance.  In  fact,  E(S^)  <  o‘.  i'his  means 
that  the  confidence  interval  generated  firom  Equation  (2)  is  narrower  than  it 
should  be;  thus  the  coverage  of  the  confidence  interval  is  much  lower  than  the 
nominal  value.  One  way  to  overcome  the  effects  of  correlated  data  is  to  batch 
the  observations. 
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The  aext  graph  demonstrates  the  effect  of  batching  the  observations  in 
successively  larger  groups  for  the  case  of  4)=0.9.  Individual  observations  are 
first  grouped  into  batches  of  m  units  each.  Figure  10  displays  batch  sizes 
where  m=2,  32,  256  observations,  the  mean  value  of  each  of  these  groups  is 
then  used  to  represent  that  group.  lOO  replications  of  sixteen  such 
observations  are  made  and  the  resulting  means  and  standard  errors  are  plotted 
for  each  value  of  m. 


Figure  10:  AR(1)  Batched  in  sizes  m=2,  m=32,  m=256 

with  16  batches 

The  graph  presented  in  Figure  lO  visually  provides  the  manager  with 
insight  into  the  effects  of  batching  correlated  data.  The  value  of  the  mean  is 
plotted  on  the  x-axis  against  the  value  of  the  resxjltant  standard  deviation  on 
the  y-axis.  One  can  see  that  the  data  when  grouped  by  2's  is  stiU  widely 
dispersed  and  coverage  is  poor.  As  the  batch  size  increases  an  obvious  pattern 
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stfiurts  to  become  clear,  coverage  improves  as  the  values  cluster  within  the 
confidence  interval  and  the  standard  error  values  decrease  as  the  sample  size 
increases.  Eventually  we  deserve  the  pattern  that  we  expect  t'rem  i.i.d.  normal 
process. 

The  decision  maker  can  use  such  graphical  presentations  to  help  him 
identify  the  trends  involved  and  use  them  to  improve  his  decision  making.  By 
batching  the  individual  data  points  the  correlation  of  the  data  diminishes  and 
is  reflected  in  the  clust(!ring  of  the  data  for  batches  of  increasing  size. 

Graphs  for  four  additional  cases  using  the  AR(1)  model  can  be  found  in 
Appendix  B.  The  graphs  present  the  effects  on  correlation  for  various  values 
of  (J).  The  cases  include  4)=0.9,  4)=0.5,  (1)=0.0.  4)=-0.5,  and  (l)=-0.9.  The  graphic 
presentation  in  each  case  is  similar  to  that  presented  above.  By  analyzing  and 
comparing  the  individual  graphs  a  manager  can  begin  to  understand  the 
significance  of  the  correlation  in  his  data. 

The  cases  where  (()=-0.5  and  (t)=-0.9  show  the  effect  of  negative  correlation 
of  data.  While  the  values  generated  in  these  cases,  even  for  small  batch  sizes, 
tend  to  provide  higher  than  nominal  coverage,  the  standard  error  is  somewhat 
smaUer  than  for  the  positively  correlated  cases,  where  (j)  is  greater  than  0. 

The  AR(1)  model  with  <})=0.6  shows  improved  coverage  at  all  batch  sizes 
when  compared  with  the  (j)=0.9  example  shown  above.  This  is  because  (J)  is 
smaller  and  consequently  each  successive  data  point  is  less  correlated  to  its 
predecessor  than  in  the  0.9  case.  The  improvements  in  coverage  afforded  by 
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batching  the  data  continue  to  apply  and  the  graphs  clearly  demonstrate  that 
point. 

The  AfUl)  model  with  4>=0.0  is,  in  fact  the  Normal  (0,1)  case.  By  reducing 
the  correlation  valu3  to  zero  there  is  no  correlation  with  the  preceding  values 
and  the  graphs  depict  the  relationship  that  would  be  expected  under  the 
assumptions  of  the  classical  model:  the  coverages  are  approximately  the  same 
as  the  nominal  value  (90%  in  this  case),  and  the  standard  error  decreases  as 

m  (batch  size)  gets  larger,  at  a  rate  of  l/v'S  . 

The  spreadsheet  designed  to  provide  the  graphical  output  in  this  section 
can  be  found  in  Appendix  B.  Macro  codmg  is  used  to  replicate  the  calculations 
of  4,096  individual  variables  100  times  in  the  generation  of  each  graph 
presented.  Each  run  consists  of  generating  the  individual  values,  grouping 
those  into  sequential  samples  of  vaunous  sizes,  and  plotting  the  mean  emd 
standard  deviation  of  sixteen  like  samples  on  a  graph.  Spreadsheet 
programming  has  thus  condensed,  millions  of  calculations  into  simple  visual 
presentations  that  can  have  a  major  impact  on  the  decisions  confronting  a 
manager  dealing  with  such  correlated  data. 

C.  EAR(1)  MODEL 

The  EAR(1)  process  is  defined  as: 
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vitA  probability:  4>  ,4) 

=  ■*■  t/ith  probability:  l-<t» 

where  -  i.i.d.  Exp(l). 

The  EAR(1)  model  has  the  same  autocorrelation  structure  as  in  the  AR(1) 
model,  except  non-normal  (exponential)  error  terms  are  added  It  is  worse 
than  the  AR(1)  model  in  terms  of  coverage  for  a  given  value  of  (J).  This  is  due 
to  asymmetric  error  terms.  The  resultant  variables  show  a  steady  pattern 


interrupted  periodically  by  sudden  and  dramatic  jumps.  This  type  of  pattern 
is  representative  of  situations  such  as  machinery  repair  in  which  the  operation 
of  a  machine  steadily  declines  until  maintenance  or  calibration  is  performed 
and  the  pattern  of  wear  or  maladjustment  starts  over  Gigain.  Figure  11 
provides  an  example  of  this  type  of  pattern  based  on  100  sequential 
observations  using  (t>=0.9. 
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To  further  explore  the  properties  of  an  EAIKl)  model  a  number  of  points 
are  again  generated  and  graphed  as  in  the  AR(1)  model  above.  Figure  12  shows 
lOO  points,  each  representing  the  sample  mean  and  sample  standard  deviation 
of  a  series  of  16  sequential  values  individually  generated  from  Equation  (4). 
This  graph  shows  that  the  values  are  widely  dispersed  and  highly  asymmetric. 
Again  the  coverage  of  the  true  mean  (p=l  in  this  case)  is  much  smaller  than 
the  nominal  value  0.90.  In  fact  only  25%  of  the  confidence  intervals  actually 
cover  the  true  mean. 


Batching  can  help  overcoming  both  the  effects  of  data  correlation  and  non¬ 
normality  inherent  in  the  EAR(1)  model.  The  next  graph  demonstrates  the 
effect  of  batching  the  values  in  successively  larger  groups  for  the  case  of  (|)=0.9. 
Individual  observations  are  first  grouped  into  16  batches  of  size  2, 32,  and  256 
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the  mean  value  of  each  of  these  groups  is  then  used  to  represent  that  group. 
100  replications  of  sixteen  such  observations  are  made  and  the  resulting  means 
and  standard  errors  are  plotted  in  Figure  13. 


0  12  3  4 

0  »25  4  %  325  «  »2BB5  - mo  I 

Figure  13:  EAR(l)  <|k=0.9  Batched  Values  m=2,  m=32,  m=256 

Again  the  decision  maker  can  easily  see  firom  the  graphical  presentation 
that  as  the  batch  size,  m,  is  increased  the  resulting  values  tend  to  cluster  more 
tightly  and  the  coverage  improves.  The  individual  points  plotted  display 
greater  symmetry  and  the  standard  errors  are  progressively  reduced  with  each 
increase  in  batch  size,  leading  eventually  to  the  acceptability  of  using  classic 
confidence  interval  procedures  to  describe  the  behavior  of  some  of  the  larger 
size  batches. 

The  spreadsheet  programs  and  additional  graphs  for  the  EAR(1)  model 
can  be  found  in  Appendix  C.  The  spreadsheet  programming  for  the  EAR(1) 


36 


model  is  virtually  identical  to  that  in  the  AR(1)  model.  According  to  Sargent. 
Kang  and  Goldsman  (1992),  the  EAR(1)  model  requires  larger  batch  size  to 
satisfy  the  assumptions  required  for  classical  confidence  interval  procedures, 
than  the  AR(1)  model  due  to  its  non- normal  error  terms.  The  additional 
graphs  provided  consider  the  cases  where  (1>=0.9,  (j)=0.5  and  (J)=0.0.  EAR(1) 
model  with  (j>=0.0  is  the  i.i.d.  Exp(l)  case. 

D.  hVHl  QUEUING  MODEL 

The  M/M/1  queuing  model  represents  a  single  channel  waiting  line  model 
with  Poisson  arrivals  and  exponential  service  times.  The  model  is  designed  to 
study  the  waiting  time  for  successive  customers  in  the  queue.  The  following 
equations  (from  Anderson,  Sweeney,  WiUiams,  1991)  specify  this  relationship: 


(4) 


p(p^-A) 


(5) 


We  used  the  values,  A,=0.8,  n=1.0  tor  this  case.  In  a  system  such  as  this  the 
values  generated  for  W,  are  obviously  dependent  on  values  immediately 
preceding.  Balking  is  not  considered  and  the  model  confines  itself  to  only  those 
customers  that  actually  join  the  queue.  A  graphic  depiction  of  the  situation  is 
presented  in  Figure  14,  where  the  horizontal  axis  represents  each  individual 
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customer  that  sequentially  enters  the  waiting  line  and  the  vertical  axis  shows 
the  customer  waiting  time. 


Figure  14:  100  Observations  from  the  M/M/1  Model 

As  shown  in  Figure  14  the  variables  generated  in  an  M/M/l  process 
appears  to  be  very  dependant  on  their  predecessors.  In  fact  such  a  system  is 
considered  to  be  both  highly  positively  correlated  and  non-normal.  Again  we 
plot  lOO  points,  each  representing  the  sample  mean  and  sample  standard  error 
of  a  series  of  16  sequential  values  individually  generated  from  Equations  (4) 
and  (5). 
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Figure  15  shows  the  values  to  be  widely  dispersed,  and  asymmetric.  The 
coverage  is  poor,  only  22  confidence  intervals  generated  fi^m  the  100 
replications  actually  cover  the  true  mean  (Wq=4  in  this  case).  Again  this  poor 
coverage  is  the  result  of  correlation,  the  waiting  time  of  any  individual 
customer  is  dependant  on  both  the  number  of  customers  in  line  ahead  of  him 
and  the  length  of  time  required  to  service  each  of  those  customers. 

Once  again  batching  can  help  overcoming  both  the  effects  of  data 
correlation  and  non-normality,  even  when  the  effects  are  quite  large  as  in  the 
NVWl  model.  The  next  graph  demonstrates  the  effect  of  batching  the  values 
in  successively  larger  groups.  Individual  observations  are  this  time  grouped 
into  only  two  batches  of  4  and  266  observations,  the  mean  value  of  each  of 
these  groups  is  then  used  to  represent  that  group.  Sixteen  such  observations 


38 


are  made  and  the  resulting  mean  and  standard  deviation  are  plotted  in  Figure 
16. 


Figure  16:  M/M/1  Batched  Values  m=4,  m=256 


Again  the  graphic  presentation  makes  it  easy  to  observe  the 
improvements  in  coverages  due  to  increasing  batch  size.  The  individual  points 
plotted  display  greater  symmetry  in  the  case  of  batch  size  256,  compared  to  a 
case  of  batch  size  4,  and  the  standard  errors  are  considerably  reduced. 

The  spreadsheet  programs  used  in  the  NVWl  model  can  be  formd  in 
Appendix  D.  Coding  for  the  M/M/1  model  is  somewhat  different  firom  the  two 
preceding  examples  and  execution  is  somewhat  slower.  The  same  spreadsheet 
and  graphical  output  formats  have  been  maintained,  and  thus  the  resulting 
output  is  easily  compared  with  that  generated  for  the  AR(1)  and  EAR(1)  cases. 
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B.  CONCLUSION 


The  three  models  discussed  in  this  chapter  are  representatives  of  a 
number  of  problems  that  the  logistics  decision  maker  may  face.  Often  these 
problems  are  treated  as  if  the  variables  involved  were  independent.  While  this 
undoubtedly  the  case  for  some  problems,  it  is  not  always  true. 

Unrecognized  correlation  of  data  can  seriously  affect  the  validity  of  a 
managers  decisions.  Such  correlation  can  not  be  eliminated,  because  it  is  a 
reflection  of  the  natural  processes  and  inherent  characteristics  of  systems.  The 
problem  facing  decision  makers  is  how  to  identify  and  deal  with  such 
correlation. 

Classical  confidence  interval  procedures  are  valid  only  if  independent, 
identically  distributed  normal  assumptions  are  satisfied  This  may  not  always 
be  the  case.  Various  logistics  problems  such  as  turn  around  time  for 
equipment  repair,  warehousing  and  transportation  problems  aU  may  be  highly 
correlated 

Spreadsheet  analysis  and  graphical  presentation  can  help  the  decision 
maker  visualize  such  systems  and  recognize  the  consequences  involved  While 
the  processes  discussed  in  this  thesis  are  all  based  on  simulated  data 
spreadsheet  programming  is  not  limited  to  such  simulations.  Any  historical 
data  or  physical  sampling  procedure  can  provide  the  raw  data  for  spreadsheet 
analysis.  When  this  data  is  available  it  can  provide  important  insight  into  the 
decision  makers  real  world  problems. 
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V.  PROCESS  CONTROL 


A.  INTRODUCTION 

Process  control  is  another  area  in  which  graphical  spreadsheet  analysis 
shows  considerable  promise.  Items  successively  produced  in  an  ongoing 
production  process  are  not  likely  to  be  totally  independent.  As  discussed  in 
Chapter  IV,  they  will  exhibit  both  random  and  systematic  variability 
dependant  on  the  characteristics  of  the  production  system  involved. 

This  systematic  variation  can  become  quite  problematic  in  the  area  of 
acceptance  sampling.  Sampling  plans  have  been  calculated,  published  and 
standardized.  Unfortunately,  these  standard  plans  generally  assume 
independent,  identically  distributed  observations  in  their  Construction- 
Application  of  such  plans  when  the  processes  involved  are  not  in  fact 
independent  can  lead  to  serious  consequences.  In  particular  it  increases  the 
risk  of  making  a  wrong  decision  based  an  ill  suited  sampling  plan. 

Type  I  error  is  considered  to  be  the  rejection  of  a  good  lot  during  the 
sampling  process  and  Type  II  error  is  the  acceptance  of  a  bad  lot.  Type  I  error 
is  considered  to  be  Producer's  Eiisk  as  the  producer  bears  the  cost  of  choosing 
to  reject  a  good  lot.  Type  II  error  is  considered  to  be  the  Consumers  Risk 
because  it  is  the  chance  that  a  bad  lot  will  not  be  detected  by  the  sampling 
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plan  employed.  Attempts  to  balance  acceptable  risk  levels  of  each  type  further 
compound  the  problem. 

In  some  situations  high  producer's  risk  is  acceptable  or  oven  desirable 
because  of  serious  repercussions  that  could  occur  hrom  passing  even  1  bad  lot. 
In  cases  such  as  weapons  manufacture  or  the  production  of  nuclear  power 
components  the  danger  of  a  component  not  working  properly  is  so  significant 
that  the  consumer's  risk  must  be  reduced  till  it  is  virtually  non-existent.  This 
is  a  costly  preposition  for  the  producer,  but  it  is  generally  considered  more 
desirable  than  the  alternative. 

The  opposite  situation  is  equally  likely  to  exist  however.  For  cheap, 
disposable,  harmless  items  the  consumer  is  often  willing  to  put  up  with  more 
risk  at  his  level.  If  a  fifty  cent  ink  pen  doesn't  work  properly  he  is  Mailing  to 
buy  another  that  may  or  may  not  work  any  better,  rather  than  buy  a  two 
hundred  dollar  Mount  Blanc  model  that  he  is  sure  will  work  as  perfectly. 

All  this  consumerism  and  logic  is  however  wasted  if  the  acceptance 
sampling  plan  in  use  does  not  properly  account  for  the  process  it  is  intended 
to  monitor. 

B.  THE  PROBLEM 

In  an  ongoing  production  process  it  is  often  very  likely  that  there  is  a 
correlation  between  successive  items  produced.  Standard  acceptance-sampling 
plans  are  usually  not  designed  to  deal  with  such  correlation  and  consequently 
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can  result  in  costly  and  time  consuming  errors.  Computing  the  probabilities 
required  to  design  sampling  plans  for  general  correlated  processes  is  often 
complex,  and  sometimes  intractable.  In  his  paper.  Nelson  (1990)  proposed  a 
methodology  for  estimating  single-sampling  attribute  plans  for  any  production 
process  that  can  be  simulated. 

In  a  typical  single  sampling  plan  of  (n.c),  n  samples  are  chosen  for 
inspection,  and  c<n  is  the  acceptance  number  such  that  the  lot  is  declared 
acceptable  if  no  more  than  c  defective  items  are  discovered  in  the  sample.  The 
values  of  n  and  c  are  chosen  to  provide  a  prespecified  producer's  risk  and 
consumer's  risk.  In  this  chapter,  we  develop  a  spreadsheet  program  to 
understand  the  pitfalls  of  using  standard  sampling  plEins  when  the  process  in 
fact  is  correlated. 

First  we  start  with  i.i.d-  measiirement  (i.e.  (j)=0.0)  for  which  standard 
acceptance  sampling  plans  are  appropriate.  We  use  process  X  as  the  in¬ 
control-process  that  has  the  normal  distribution  with  a  mean  of  10  and  the 
variance  1.  The  process  X'  is  considered  an  out-of-control  process  that  has  the 
same  variance  as  X,  but  a  process  mean  that  has  shifted  to  11.294. 

For  such  processes  a  defective  item  is  defined  as  an  item  from  the 
process  that  is  beyond  the  tolerance  limits.  In  this  example  we  chose  the 
upper  tolerance  limit  Tu=12.5768  and  the  lower  tolerance  limit  Tl=7.4242  as 
shown  in  Table  1.  These  numbers  are  chosen  so  that  the  producer's  risk  is 
O.Ol  and  the  consumer's  risk  is  0.10;  i.e. 
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Pr  (li  <  <  Xy)  =0.99  and 


Pt  (x^  <  Jf'  <  Xp)  =0.10. 


1  quality 

4 

a 

Tl 

P 

AQL 

10 

1 

7.4242 

12.5758 

0.01 

LTPD 

11.2940 

1 

7.4242 

12.5758 

0.10 

rarr 


AQL  in  Table  I  stands  for  Acceptable  Quality  Level  and  is  the  desired 
maximiim  probability  of  stopping  the  production  process  if  it  is  operating  at  the 
acceptable  level.  LTPD  stands  for  Lot  Tolerance  Percent  Defective  and  is  the 
minimum  probability  of  rejecting  the  process  when  it  is  operating  at  the 
unacceptable  level.  Now  using  n=39  and  c=l,  we  generate  39  samples  from  the 
process  for  X  and  X'  (Figure  17). 

From  Figure  17  we  can  see  that  4  values  exceed  the  tolerance  level  for  the 
out-of-control  process  CX*)  and  that  all  of  the  values  fall  within  the  tolerance 
limits  for  the  in-control-process.  The  graph  shows  both  processes  together  and 
can  help  the  decision  maker  visualize  the  differences  between  them. 
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Figure  18  shows  these  results,  for  a  different  replication,  in  tabular 
format.  In  this  figure  C-Lot  stand  for  the  lot  produced  by  the  in-control - 
process  and  OOC  Lot  stands  for  the  lot  produced  by  the  out-of-control  process. 
The  Defects  column  shows  how  many  items  exceeded  the  tolerance  limits  for 
each  respective  lot  and  the  Lot  column  provides  a  recommendation  based  on 
the  established  decision  rule,  reject  if  c>l. 


Defects 

Lot 

Cixi 

1 

Axsct 

OOC  Lot 

3 

Reject 

Figure  18:  Accept  /  Reject  Table 


We  repeated  this  process  1,000  times  and  obtEuned  the  results  shown  in 
Figure  19. 
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C-Lot 

OOC_ljot 


Total  Accept  Total  Reject 
948  52 

81  919 

Figure  19:  Result  of  1000  Trials 

Expected  values  for  this  process  can  be  calculated  from  the  equations: 

Pzoducei'a  risk  =  l  -  =  0.0581 

jc-oV  .  p  ^ 

Conaumei'a  riak  =  0 . 10-^  0 . =  0.0876 

Table  II  shows  the  effect  of  increased  correlation  as  4>  is  systematically 
increased. 


u 

Producer's  Risk 

Consumer's  Risk 

1  Theoretical  Values 

.0581 

.0876 

0.0 

.062 

.081  1 

0.1 

.059 

.110 

0.3 

.064 

.128 

.06 

.072 

.170 

0.7. 

.087 

.267 

0.9 

.082 

.496 

Table  II. 
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For  the  4)=0  ('ase,  the  computer  results  are  close  to  the  theoretical  values, 
as  they  should  be.  However,  as  you  change  the  (j>  values  to  4)=0.1,  0.3,  0.6,  0.7, 
0.9  you  increase  the  correlation  in  the  process  while  keeping  the  process 
variance  as  is  (i.e.  a"=l),  both  the  producer's  risk  and  the  consumer's  risk 
increasingly  differ  from  the  i.i.d  case.  As  seen  in  Table  II  the  correlation  in  the 
process  causes  a  problem,  the  risk  estimators  are  significantly  different  from 
the  theoretical  expectation  of  the  normal  i.i.d.  case. 

C.  CONCLUSION 

This  chapter  has  demonstrated  how  correlation  can  have  a  negative 
impact  on  sampling  plans.  The  decision  maker  needs  a  tool  to  estimate  the 
risk  associated  with  such  plans.  Spreadsheets  are  an  excellent  tool  that  can 
provide  the  calculation  power  necessary  and  the  graphical  ability  to  help  the 
decision  maker  better  understand  the  problem.  Complex  relationships  can  be 
clarified  and  the  manager  is  provided  with  graphs  and  tables  to  help  him  to 
draw  conclusions  that  may  not  be  otherwise  obvious. 
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VI  CONCLUDING  REMARKS 


This  thesis  shows  through  the  use  of  several  diverse  logistics  problems 
that  spreadsheets  can  provide  a  good  foundation  for  logistics  decision  support 
systems.  The  availability  of  the  personal  computer  and  spreadsheet  programs 
make  them  a  logical  choice  for  such  systems.  Further,  the  simplicity  of  their 
operation  and  the  powerful  output  options  available  add  greatly  to  their  value 
for  this  pxirpose. 

We  have  shown  that  spreadsheet  models  can  be  developed  to  provide 
decision  support  for  a  broad  range  of  logistics  problems.  Models  were 
presented  in  the  areas  of  inventory  simulation,  time  series  analysis,  data 
correlation  and  process  control,  and  the  value  of  graphical  presentation  as  an 
aid  to  decision  making  for  each  case  was  discussed.  Finally  we  have  shown 
that  spreadsheet  programs  can  provide  the  logistics  decision  maker  with  the 
tools  he  needs  to  develop  programs  himself  to  solve  the  daily  problems  which 
confront  him.  The  value  of  spreadsheet  analysis  as  a  foundation  for 
logistics  decision  support  lies  more  in  its  ready  availability  and  simplicity  than 
in  its  computational  power.  Custom-made  programs,  professionally  developed 
to  handle  specific  problems  would  likely  be  far  more  capable  and  could  provide 
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far  more  elaborate  graphical  presentations.  While  this  might  be  preferable 
from  a  scientific  point  of  view  it  would  be  difficult  to  implement. 

The  cost  and  difficulty  involved  in  the  development  of  specialty  programs 
would  severely  limit  the  applications  which  could  be  automated.  The 
logistician  needs  systems  that  allow  him  to  quickly  model  a  system  in  which 
he  is  interested  and  obtain  timely  results  to  help  him  make  decisions.  While 
spreadsheets  may  not  provide  the  most  elegant  decision  support  system  they 
certainly  provide  sufficient  power  and  flexibility  support  many  of  the  decisions 
that  confront  a  logistician  on  a  daily  basis. 

The  graphical  options  available  in  spreadsheet  packages  offer  the  decision 
maker  with  a  broad  array  of  display  options.  These  presentations,  previously 
costly,  time  consuming  and  difficult  to  obtain  are  now  available  to  the  user 
with  a  few  keystrokes.  Spreadsheet  programming  offers  the  logistician  a 
desktop  solution  to  the  complex  logistics  problems  that  confront  him  on  a  daily 
basis. 
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APPENDIX  A  INVENTORY  MODEL  CODE 


A.  MAIN  MACRO 
ftmctioii  NewValues 

WelcomeAns  =DIALOG.BOX(dWelcomeScreen) 

=  IF(WelcomeAns=8,GOTO(B5),) 

=  IF(WelcomeAiis=FALSE,CLOSE(FALSE)) 

InputAns  =DIALC)G.BOX(dInputSheet) 

=  IF(InputAns=FALSE,GOTO(El)) 

=OPEN("mam.xls") 

=OPEN("chartLxlc") 

=RUN(Initialization) 

=FORMULA(DAY(NOW0)&"  /  ''&MONTH(NOW0)&"  / 
”&YEAR(NOW0)&'’("&HOUR(NOW0)&":'’&MINUTE 

(NOW0)&")",R40) 

=RUN(Headers) 

=RUN(Repeater) 

=RUN(PlotNetInv) 

=RUN(PlotNetInvUpdate) 

Final  Ans  =DIALOG.BOX(dFinalAns) 

=RUN(FinalScreenAns) 

=RETURN0 


ftmctxm  FinalScreenAxis 


=IF(FinalAns=l) 

=  ACTIVABi:"chartl.xlc") 

=  UNHIDEC'chartLxlc") 

Keep  changmg=  ON.KEY("~",’'work.xlm!R47c2") 
=  PAUSED 
=  ON.KEYC'~",) 

=  HIDEO 
=  GOTO(Bl5) 
=ELSE.IF(FinalAns=lO) 
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=ELSE.IF(FitialAns=l  1) 

=  ACTIVATE("mam.xls") 

=  UNHIDE(''main.xJa“) 

=  ON.KEY(”~''."work.xlm!R47c2'') 
=  PAUSEO 
=  ON.KEYC'-",) 

=  HIDEO 

=  GOTCXBIS) 

=ELSE.IF(FiiialAns=l2) 

=ELSE.IF(FiiiaIAns=l3) 

=ELSE.IF(FinalAns=FALiSE) 

=RUN(Clo8eAllWorkFaes) 

=END.IF0 

=RETURNO 


ftmction  MacGofntmue 

=RESUME(1) 

=RETURN0 


function  aoaeAUWorkFiies 

=ACTI  VATEC'chartl  .xlc") 

=SAVE.AS("D:\THESIS\invmod\LASTWORK.XLC". 

l."",FALSE."".FALSE) 

=ACTIVATE("MAIN.XLS") 

=SAVE.AS("D:\THESIS\mvmod\LASTWORK.Xls", 

1,"", FALSE,"", FALSE) 

=ACTlVATE("LASTWOElILXLS") 

=SELECTC'Rlcl:R"&WORK.XLM!Repetition&"c2l".  "RlCl") 
=SELECT("R1:R16384") 

=CLEAR(3) 

=SELECT("R1C1") 

=SAVE.AS("D:\THESIS\invmod\MAIN.XLS",l,"",  FALSE,"", 

FALSE, 

=CLOSE0 

=ACTIVATE("LASTWORK.Xlc") 

=SAVE.ASrD:\THESIS\mvmod\chartl.xlc", 

1,"",FALSE,"",FA1SE) 

=CLOSE(TRUE) 

=RETURN0 
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ftmction  Initialization 


=ACTIVATE("MAIN.X1JS") 

=SELECT("mairLxls!Rlcl;R"&WORK.XLM! 

Repetition&"c2 1  ","Rl  C 1  ”) 

=CLEAR(3) 

=ACTIVATE("WORK.XLM") 

=FORMULA(0,X3) 

=FORMULA(l  .Repetition) 

=FORMULA(0,AC1) 

=FORMULA("".X5:AA7) 

=FORMULA("‘'.X9) 

=SELECT(WORK.XLM!  AC3:  AE 1 96) 

=  CLEAR(3) 

=FORMlJLA(99999999999,X7) 

=FORMULA(99999999999.X6) 

=SELECT(WORK.XLM!AC3:AE398) 

=  CLEAR(3) 

=FORMULA(0,AA5) 

=FORMULA(0,AA6) 

=FORMULA(0,AA7) 

=RETURN0 


function  Headers 

=FORMULA("Event'’;D:\THESIS\PARTl\INVMOD\MAIN. 

XI^’!K2) 

=FORMULA("Event'’;D:\THESIS\PARTl\INVMOD\MAIN. 

XLS'lLl) 

=FORMUlJV(''Time", 'D:\THESIS\PARTl\INVMOD\MAIN. 
XLS'!L2) 

=FORMULAC'Sys'', 'D:\THESIS\PART1\INVM0D\MAIN. 
XLS'!M1) 

=FORMULA("Recpt", 'D:\THESIS\PARTl\INVMOD\MAIN. 
XLS’!M2) 

=FORMULAC'Cust";D:\THESIS\PARTl\INVMOD\MAIN. 

XLS’INl) 

=FORMULA("Req'', 'D:\THESIS\PARTl\INVMOD\MAIN. 
XLS'!N2) 
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=FORMULA("Sys",'D:\THESIS\PARTl\INVMOD\MAIN. 

XLS’!01) 

=FORMULA("Order",'D:  \THES  IS \PART1  \IN  VMODNMAIN. 
XLS‘!02) 

=FORMULA("Qty'‘.'D:\THESIS\PARTl\INVMOD\MAIN. 

XLS'IPl) 

=FORMULA('’  Iss’’,’D;  \THES  IS  \PARTl  \  INVMOD  \MALV. 
XLS'!P2) 

=PORMULA(''Qty", 'D:\THES1S\PART1\INVM0D\MA1N. 
XLS'lQl) 

=FORMULAC'BO'', 'D:\THESIS\PARTl\INVMOD\MAIN. 
XLS'!Q2) 

=FORMULA(''Qty''.'D:\THESIS  \PART1  \INVMOD\MAIN. 
XLS'lRl) 

=FORMULA("OH".'D:\THESIS\PARTl\INVMOD\MAIN. 

XLS’!R2) 

=FORMULA("  Qty "  ,'D:  \THESIS  \PART1  \  IN  VMOD\M  AI N . 
XLS'lSl) 

=F0RMULA("00";D:\THESIS\PART1\INVM0D\MAIN. 

XLS'!S2) 

=FORMULA("Net",'D:\THESIS\PARTl\INVMOD\MAIN. 

XLS'lTl) 

=FORMULA(''Inv'', 'D:\THESIS\PART1\INVM0D\MAIN. 
XLS'!T2) 

=FORMULA("Inv".'D:\THESIS\PARTl \INVMOD\  MAIN. 
XLS'IUl) 

=FORMULA("Pos'' .'D:\THESIS  \PART1  \INVMOD\MAIN. 

XLS’!U2) 

=FORMULA("Beg 

Bal", 'D:\THESIS\PART1\INVM0D\MAIN.  XLS'!K4) 

=FORMULA(VALUE(WORK.XLM!OnHand);D:\THESIS\PARTl\ 

INVMOD\MAIN.XLS’!R4) 

=FORMULA(VALUE(WORK.XLM!OnOrder);D;\THESIS\ 
PARTI  \INVMOD\MAIN.XLS'!S4) 


=FORMULA(VALUE(WORK.XLM!OnBackorder),'D:\THESIS\ 

PARTI  \INVMOD\MAIN.XLS'!Q4) 
=FORMULA('D:\THESIS\PARTl\INVMOD\MAIN.XLS'!R4 
'D:\THESIS\PART1  \INVM0D\MAIN.XLS'!Q4;D:\ 
THESIS\PARTl\INVMOD\MAIN.XLS'  !T4) 
=F0RMULACD:\THESIS\PART1\INVM0D\MAIN.XLS'!R4 
'D:\THESIS\PART1  \INVM0D\MAIN.XLS'!Q4+'D: 
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\THESIS\PARTl\INVMOD\MAlN.XLS'!S4,'D:VniESlS\ 
PARTI  \INVM0D\MAIN.XIJS'!U4) 
=SELECT(WORK.XLM!Vl) 

=RETURNO 


AjDCtMm  CalcCustRec 

=IF(DemandDist=  1 ) 

=  FORMULA(0,V1) 

=ELSE.  IF(DemandDist  =2) 

=  FORMUlJ^(VLOOKUP(RAND0,TEXTREF('’r3c36:r'’& 
(AG2+4)&‘'c37’'),2),Vl) 

=ELSE.IF(DemanclDist=3) 

=  FORMULA(VLOOKUP(RAND0.TEXTREF('’r3c36;r'’& 
(AG2+4)&"c37").2),Vl) 

=ELiSE.  IF(DemandDist=4) 
minmDD  =  DemandMean 

maxmDD  =  DemandStdDev 

=  F0RMULA("'’,B146) 

=  FOR(''count'\l,ABS(miomDD-maxniDD),l) 

=  FORMULA(Bl46&"+randO",Bl46) 

=  NEXTO 

=IF(randserie8DD=’"',FORMULA('’+randO'’,Bl46)) 

=  F0RMULA(BIGH1XB146,(LEN(B146)-1)),B146) 

=  FORMULA("=''&TEXT(ininmDD,0)&'’+  C’&randseriesDD 
&T,B144) 

UDAnsDD  =0+(RAND0) 

=  FORMULA(UDAiisDD,Vl) 
randseriesDDrandO 

=ELSB.IF(DemandDist>4) 

=  FORMULA(RANDO,V1) 

=END.IF0 

=IF(RecQtyDist=l) 

=  FORMULA(0,V2) 

=ELiSE.IF(RecQtyDist=2) 

=  FORMUIJ^(VLOOKUP(RAND0,TEXTREF(''r3c33:r"& 
(AG2+4)&"c34"),2),V2) 

=ELSE.IF(RecQtyDi8t=3) 

=  FOElMUIJX(VLCXDKUP(RAND0,TEXTREF("r3c33:r''& 
(AG2+4)&''c34'’),2),V2) 

=ELSE.IF(R^QtyDi8t=4) 
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minmQD  =RecQtyMean 
maxmQD  =RecQtyStdDev 

=FORMULA("".Bl69) 

=FOR("  count" ,  1 ,  ABS(minmQD-maxmQD),  1 ) 
=FORMULA(B  169&"+randO".B  169) 

=NEXTO 

=IF(randseriesDD="",FORMULA("  +randO",B  1 69)) 
=FORMULA(RIGHT(B  1 69.(LEN(B  1 69)- 1  )),B  1 69) 
=FORMULA("  ="&TEXT(mininQD,0)&"  +("&arandseriesQD 
&")",B167) 

UDAnsQD  =2+(RAMX)+RAND0+RAND0+RAND0) 
=FORMULA(ROUND(UDAnsQD,0),V2) 
randsoriesQDrandO+randO+randO+randO 
= E  LS  E .  IF(RecQty  Dist>4) 

=  FORMIJLA(RAND0,V2) 

=END.IF0 

=FORMULA(V  1  +X3^5) 

=FORMULA(V2,Y5) 

=FORMULA("  ",Z5) 

=FORMULA(l,AA5) 

=RETURNO 


ftmction  CalcSysReoider 

=IF(GET.CEm6,TEXTREF(ADDRESS(Repetition 
+3,2 1  „FALSE,"mam.xl8")))<=TVp©InvR) 

=  FORMULA(X3.X6) 

=  FORMULA("",Y6) 

=  FORMULA(TypeInvQ+(TypeInvR-GET.CELL(5, 
TEXTREF(ADDRESS(Repetition+3,2 1  „FALSE, 
"main.xls")))),Z6) 

=  RUNCCalcPLT) 

=  FOElMULA((PLTAnswer+X6),TEXTREF(ADDRESS 
(ACl  +3,29,1,TRUE,),TRUE)) 

=  F0RMULA(Y6,TEXTREF(  ADDRESS(AC  1  +3, 

30.1, TRUE,),TRUE)) 

=  FORMULA(Z6,TEXTREF(ADDRESS(ACl+3, 

31.1,  TRUE,), TRUE)) 

=  FORMULA(ACl+l,ACl) 

=FORMULA(l,AA6) 

=ELSE0 

=  FORMULA(99999999999,X6) 
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=  F0RMULA("’',Y6) 

=  FORMULA("’',Z6) 

=  FORMULA(O.AA6) 

=END.IFO 

=RETURNO 


Buhroutine  CalcFLT 

=IF(PLTDist=l) 

=ELSE.IF(PLTDist=2) 

=ELSE.IF(PLTDi8t=3) 

=ELSE.IF(PLTDist=4) 

=ELSE.IF(PLTDist=6) 

=  FORMULA(PLTMeaii,B2l6) 
=ELSE.IF(PLTDi8t>6) 
=END.IF0 
PLT  Answer  5 

=RETURN(PLTAnswer) 


ftmction  CalcSysRecpt 

=GET.CELL(6,AC1) 

=IF(AC1=0) 

=  FORMULA(99999999999,X7) 

=  FORMULAC'.Y?) 

=  FORMULA("",Z7) 

=  FORMULA(0,AA7) 

=ELSE0 

=FORMULA.GOTO(''work.xlin!R3C3l") 

=SE1JECTC'R3C29:R20C31") 

=SET.DATABASEO 

=SORT(l,"R3C29M) 

=  IF(X7>=9999999999,C0PY(AC3:AE3,X7),G0T0(B241)) 
=  F0RMULA(1,AA7) 

=  F0RMULA('"’,AC3) 

=  F0RMULAC"',AD3) 

=  FOElMULAC"',AE3) 

=  SELECTC'R3C29:R2(X)C3l") 

=  SET.DATABASEO 
=  SORT(i;'r3c29M) 

=  FORMULA(ACl-l,ACl) 
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function 


=END.IFO 

=RETURNO 


Repeater 

=IF(AA5=0,RUN(CalcCustRec),) 

=IF(AA6=0,RUN(CalcSysReorder),) 

=IF(AA7=0,RUN(CalcSysRecpt),) 

=FORMULA(MIN(X5:X7),X9) 

=IF(X9=X5) 

=  FORMULA(W5.TEXTREF(ADDRESS(X4+4. 1 1 . 
„''main.xls‘'),TRUE)) 

=  FORMULA(X5,TEXTREF(ADDRESS(X4 
+4, 12,„"mainuxls"),TRUE)) 

=  FORMULA(Y5,TEXTREF(ADDRESS(X4+ 

4, 1 4.„”maiTLxls"),TRUE)) 
=RUN(CalcCustRecSub) 

=  FORMULA(X3+VlJ(3) 

~  FORMULA(Repetition+ 1  ,Repetitioii) 

=  FORMULAC  ".X6) 

=  FORMULAC  ".Y5) 

=  FORMULAC  ",Z5) 

=  FORMULA(O.AA6) 

=ELSE.IF(X9=X6) 

=  FORMULA(W6,TEXTREF(ADDRESS(X4 
+4,1 1  ,„"main.xls"),TRUE)) 

=  FORMULA(VALUE(X6),TEXTREF(ADDRESS(X4 
+4.  l2,„"main.xIs"),TRUE)) 

=  FORMULA(VALUE(Z6),TEXTREF(ADDRESS(X4 
+4,16,„"main.xls"),TRUE)) 

=  RUN(CalcSysfleordSub) 

=  FORMULA(X6,X3) 

=  FORMULA(Repetition+l, Repetition) 

=  FORMULA(99999999999,X6) 

=  FORMULAC '‘,Y6) 

=  FORMULAC '',Z6) 

=  FORMULA(0,AA6) 

=ELSE.IF(X9=X7) 

=  FOElMULA(W7,TEXTREF{ADDRESS(X4+4, 

1 1  ,„"mam.xls"),TRUE)) 

=  F0RMULA(X7,TEXTREF(ADDRESS(X4+4,12 
,„''main.xls"),TRUE)) 
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=  F0RMULA(Z7.TEXTREF(ADDRESS(X4+4.13 
,„"inam.xls"),TRUE)) 

=  RUN(CaIcSysRecptSub) 

=  FORMULA(X7.X3) 

=  FORMULACRepetition+l.Rejj^tition) 

=  FORMULA(99999999999.X7) 

=  FOEIMULAC  ”,¥7) 

=  FORMULAC  “,Z7) 

=  FORMULA(0,AA7) 

=END.IF0 

=IF(X3<EndTime,GOTO(B246),) 

=RETURNO 


subroutine  CalcSysReordSub 

=FORMULA(GET.CELL(5,TEXTREF(  ADDRESS 
(Repetition+3, 1 7„FALSE,"inaiiLxls"))), 
TEXTREF(ADDRESS(Repetition+4, 1 7„FALSE. 
"main-xls"))) 

=FORMULA(GET.CELL(5,TEXTREF(ADDRESS 
(Repetition+3, l8„FALSE,"main.xIs’’))),TEXTREF 
(ADDRESS(Repetition+4,18„FALSE,"mam.xls"))) 

=FORMULA(GET.CELL(5,TEXTREF(ADDRESS 
(Repetition+3, l9„FALSE,"mam.xls")))+GET. 
CELL(6,TEXTREF(ADDRESS(Repetition+4,l5. 
.FALSE, "main.xl8'’))),TEXTREF(ADDRESS( 
Repetition+4, 1 9„FALSE,''main-xls"))) 
=FORMULA(("  ="  &ADDRESS(Repetition+  4,18, 

, FALSE, "mam.xls'')&""&ADDRESS(Repetition 
+4,l7„FALSE,"main.xl8")),TEXTREF(ADDRESS 
(Repetition+4,20.,FALSE,"main.xls"))) 
=FORMULA(("="&ADDRESS(Repetition+4,l8, 

,FALSE,"main-xls”)&""&ADDRESS(Repetition 
+4,l7„FALSE,''maiiLxl8")&"+”&ADDRESS( 
Repetition+4, 1 9.,FALSE,"main.xls")),TEXTOEF 
(ADDRESS(Repetition+4,2l„FALSE,"main.xl8"))) 
=RETURN0 


suhrofutine  CalcChistRecSub 


58 


=IF(GET.CELU5,TEXTREF(ADDElESS(Repetition+3. 

l8„FALSE.''main.xls")))>=GET.CELL 

(5.TEXTREF(ADDRESS(Repetition+4, 

14.,FALSE,"main.xJs")))) 

=  FORMULA(GET.CELU5.TEXTREF(ADDRESS 
(Repetition+4,14„FALSE."main.xls"))), 
TEXTREF(ADDRESS(Repetition+4, 1 6. 

.FALSE, "main.xls"))) 

=  F0RMULA(GET.CELL(5.TEXTREF(ADDRESS 
(Repetitioa+3.l7.,FALSE,’'main.3ds"))), 
TEXTREF(ADDRESS(Repetition+4,l7, 
,FALSE,"main.xls"))) 

=  FORMULA(GET.CEm5.TEXTREF(ADDRESS 
(Repetition+3,l8..FALSE.’’inain.xJs"))). 
GET.CELU6,TEXTREF(  ADDRESSCRepet  ition+4 .14. 
,FALSE."inain.xls’'))).TEXTREF(ADDRESS 
(Repetition+4,l8„FALSE,"main.xls''))) 

=  F0RMULA(GET.CELL(5.TEXTREF(ADDRESS 
(Repetition+3.19„FALSE,''mairLxls"))), 
TEXTREF(ADDRESS(Repetition+4, 1 9,.FALSE. 
"maiiLxls"))) 

=ELSEO 

=  F0RMULA(GET.CELL(5.TEXTREF(ADDRRSS 
(Repetition+3, 1 8„FALSE,"main.xls"))), 
TEXTREF(ADDRESS(Repetition+4.l6., FALSE. 

"inaiii.xls"))) 

=  F0EiMULA(GET.CELL(6.TEXTREF(ADDRESS 
(Bepetition+3, 1 7„FALSE."main.xls"))) 
+GET.CELL(5,TEXTREF(ADDRESS(Repetition+4.14. 
.FALSE."main.xls")))-GET.CELL(5, 

TEXTREF(  ADDRESS(Repetition+3, 1 8, 

, FALSE, "maiiLxls"))),TEXTREF(  ADDRESS 
(Repetition+4, 1 7„FALSE,"inain.xls"))) 

=  FORMULA(0,TEX'niEF(ADDRESS(Repetition+4, 18. 
,FALSE.''inam.xls"))) 

=FORMULA(GET.CELL(6,TEXTREF(ADDRESS(Repetition+3, 
19„FALSE,''main.xls"))),TEXTREF(ADDRESS 
(Repetition+4, 1 9„FALSE,"m8diLxls''))) 

=END.IFO 

=FORMlJLA(("="&ADDRESS(Repetition+4, 

18„FALSE,"main.xIs'')&'’”&ADDRESS 

(Repetition+4,l7„FALSE,"main-xls")), 
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TEXT.^EF(ADDRESe(Repetition+4.20,.FALSE. 

"main.xls"))) 

=FORMULA(("=‘‘&ADDRESS(RepctiLion-r4.l8. 
.false;  main.xls")&""&ADDRESS 
(Ropetition+4,l7.,FALSE."main.)ds'’)& 

" +"&ADDRESS(Repetition+4 .19, 

, FALSE, "main.xls'')),TEXTREF(ADDRESS 
(Repetition+4,2l..FALSE."main*xls"))) 
=RETURNO 


subroutine  CalcSysRecptSub 

=IF(GET.CEIJa6.TEXTREF(ADDRESS(Ropetition+4. 

1  FALSE."main.xls")))>=CET.CELL(5. 

.  •aiiEF(AnDRESS(Rept4ition+3.l7..FALSR. 
"inam.xls")))) 

=  FOElMULA(0,TEXTREF(ADDRESS(Repetition+4. 
l7..FAlJ5E."maiiLxls"))) 

=  FORMULA(GET.CELL(5,TEXTREF(ADDRESS( 
Repetition+3,l8„FMJSE,"main.xls")))+GET. 
CELL(5,TEXTREF(ADDRESS(Repetition+4, 
l3„FALSE,"maip.xIs")))GET.CELl/5,TEXTREF 
(ADDRESS(Rep^tition+3.l7..FALSE, 
"mam.xls"))),TEXTREF(ADDRESS(Repetition+4, 
l8„FALSE.”main.xls"))) 

=  FORMULA(GET.CELL(5,TEXTREF(ADDRESS 
(Rep'>tition+3,l9„FALSE,'’main.xls''))K»ET.CELL 
(5,TEXTREF(ADDRESS(Repetition+4,l3, 

, FALSE, "mam.xl8"))),TEXTREF(ADDFESS 
(RQpetition+4, 1 9„F  ALSE,"  m  ain-xls"))^ 

=EL3E0 

=  FORMULA(GET.CELL(6,TEXTREF(ADDRESS 
(Repetition+3,l7,.FALSE,"mam-xls'')))  - 
GET.CELL(6,TEXTREF(ADDRESS(Repetition 
+4,l3.,FALSE,"main.xl8''))).TEXTREF( 
ADDRESS(Repetition+4,l7„FALSE,"inain.xls"))) 

=  FORMULA(GET.CELL(6,TEXTREF(ADDRESS(Repetition 
+3.l8„FALSE,"main.x]s"))),TEXTREF(ADDRESS 
(Repetition+4,l8„FALSE,''main.xl8'‘))) 

=  FORMLJLA(GET.CELL(5,TEXTREF(ADDRESS(Repetition 
+3,l9„FALSE,"main-xls"))X:iET.CELL(5, 
TEX'niEF(ADDRESS(Repetition+4,13.,FALSE, 
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"mairLxls''))).TEXTREF(ADDRESS(Re  petition 
+4,l9.,FALSE."main.xls"))) 

=EiND.IFO 

=FOElMULA((’'=‘'&ADDRESS(Repetition+4.18..FALSE, 
"main.xls")&''"&ADDRESS(Repetition+4. 1 7. 
,FALSE,"main.xls“)),TEX'rREF(ADDRESS 
(Repetition+4,20„FALSE,"mairLxIs"))) 
=FORMULA(("="&ADDRESS(Repetition+4.l8., FALSE. 
"main.xls")&"-"&ADDRESS(Repetition+4,l7, 

,  F  ALS  E , "  main.xls"  )&" + "  &  ADDRE  S  SCRepetition 
+4,l9„FALSE,"main.xls")).TEXTREF(ADDRESS 
(Repetition+4,2 1  „FALSE,"main.xis"))) 
=RETURNO 


function  FlotNetlnv 
nmi 

=''main.xls!r4cl2;R"&(WOElK.XLM!Repetition 
+4)&"cl2,main.xls!r4c20:r’'&(WORK.XLM! 
Repetition+4)&"  c20" 

=ACTIVATE("MAIN.XLS") 

=COPY0 

=CREATE.0BJECT(6,"R1C1",21.7.5.'’R15C9", 27.75, 

6.76,1, TRUE) 

=CHART.WIZARD(TRUE,TEXTREF(WORK.XLM!nmi),3,8,2, 
1,2,2, "Net  Inventory","Time","Units","") 
=CHART.WTZARD(TRUE,"MAIN.XLS!R4C12:R100C12. 
R4C20:R100C20"  .3,8,2, 1 ,2,2."Net  Inventory" , 
"Time"."Units","") 

=RETURN0 


function  FlotNetlnvUpdate 

=UNHIDE("MAIN.XLS  Chart  1") 

=WINDOW.MAXIMIZE0 

PlotNetInvX="niain.xls!r4cl2:r"&(Repetition+4)&"cl2" 
PlotNetInvY="main-xls!r4c20:r"&(Repetition+4)&"c20" 
=ACTIVATE("MAIN.XLS  Chart  1") 
=EDIT.SERIES(l,"NetInv",TEXTREF(WORK.XLM! 

PlotNetInvX),TEXTREF(WORK.XLM!PlotNetInvY)  „  1) 
=RETURN0 
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function  MakeDistributionTables 


=SELECT(WORK.XLM!AG;AK) 

=CLEAR(3) 

=SELECT(WORK.XLM!Vl) 

=IF(WORK.XLM!RecQtyDist=2.RUN(MakeDmdQty 

TabNormal),) 

=IF(WORK.XLM!RecQtyDist=3,RUN(MakeDmdQty 

TabPoisson),) 

=IF(WOElK.XLM!DemandDist=2,RUN(MakeDmdFreq 

TabNormal),) 

=IF(WORK.XLM!DemandDist=3.RUN(MakeDmdFreq 

TabPoisson),) 

=RETURN0 


ftmction  MakeDmdQtyrabNormal 

=FORMULA(l,AGl) 

=IF(T(TEXTREF(’'r"&AGl+ l&"c33"))=""  ,GOTO(B366), 
GOTO(B370)) 

=  FORMULAC  ","r"&AGl+l&"c33") 

=  FORMULAC  ",''r"&AGl+l&"c34") 

=  FORMULA(AGl+l,AGl) 

=  GOTO(B366) 

=FORMULA("Normal  Qty'’,AGl) 

=FORMULA(0,AG2) 

=FORMULA(0,AG3) 

=FORMULA(0,AH3) 
=FORMULA(AG2,''r’'&AG2+4&"c34") 
=FORMULA0^ORMDIST(AG2,RecQtyMe£in,RecQty 
StdDev,TRUE),"r''&AG2+4&"c33") 
=IF(VALUE(TEXTREF("r"&AG2+4&"c33"))=  1  ,GOTO 
(WORK.XLM!B379),GOTO(WORK.XLM!B377)) 
=FORMULA(AG2+ 1 ,  AG2) 

=GOTO(B374) 

=RETURN0 


ftmction  MakeDmdQtyTabPoisson 
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=F0RMULA(1,AG1) 

=IF(T(TEXTREF(’'r''&AGl  +  l&’'c33'’))=''’',GOTO(B385) 
,GOTO(B389)) 

=  FORMULAC  ",''r‘’&AGl  +  l&"c33") 

=  FORMULAC  "."r"&AGl+l&"c34") 

=  FORMULA(AG1  +  1,AG20) 

=  G0T0(B384) 

=FORMULA("Poisson  Qty".AGl) 

=FORMULA(0,AG2) 

=FORMULA(O.AG3) 

=FORMULA(O.AH3) 

=FORMULA(AG2."r"&AG2+4&'’c34") 

=FORMULA(POISSON(AG2,RecQtyMean,TRUE),"r’' 

&AG2+4&"c33'') 

=IF(VALUE(TEXTREF("r"&AG2+4&V33‘’))>=0.999999999. 

GOTO(WORK.XLM!B398).GOTO(WORK.XLM!B396)) 
=FORMULA(AG2+ 1  ,AG2) 

=GOTO(B393) 

=RETURN0 


function  MakeDmdFreqTabNormBl 

=FORMULA(l,AJl) 

=IF(T(TEXTREF("r"&AJ  1  +  l&"c36’’))="''.GOTO 
(B404),GOTO(B408)) 

=  FORMULAE  "."r"&AJl+l&"c36") 

=  FORMULAC’  ","r"&AJl+l&"c37’’) 

=  F0RMULA(AJ1+1.AJ1) 

=  GOTO(B403) 

=FORMULA("Normal  Freq”,AJl) 

=FORMULA(0,AJ2) 

=FORMULA(0,AJ3) 

=FORMULA(0,AK3) 
=FORMULA(AJ2,"r"&AJ2+4&'’c37") 
=FORMULA(NORMDIST(AJ2,DemandMean,Deinand 
StdDev,TRUE),"r"&AJ2+4&"c36") 
=IF(VALUE(TEXTREF("r''&AJ2+4&''c36"))=  I  .GOTO 
(WORK.XLMIB4 1 7),GOTO(WORK.XLM!B4 1 5)) 
=FORMULA(AJ2+l,AJ2) 

=GOTO(B412) 

=RETURN0 
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function  MakeDmdFreqTabPoisson 


=FORMULA(l,AJl) 

=IF(T(TEXTREF("r"&/\J  1  +  ''.GOTO 

(B423),GOTO(B427)) 

=  FORMULAC  ",''r"&AJl+l&"c36'*) 

=  FORMULAC  ","r''&AJl+l&"c37") 

=  FORMULA(AJl+l,AJl) 

=  GOTO(B422) 

=FORMULA(" Poisson  Preq'',AJl) 

=FORMULA(0,AJ2) 

=FORMULA(0,AJ3) 

=FORMULA(0,AK3) 

=FORMULA(AJ2,“r"&AJ2+4&”c37”) 

=FORMULA(POISSON(AJ2.DemandMean,TRUE)."r"& 

AJ2+4&V36'') 

=IF(VALUE(TEXTREF("r"&AJ2+4&''c36'’))=l,GOTO 
(WORK.XLM!B436),GOTO(WORK.XLM!B434)) 
=FORMULA(AJ2+ 1,AJ2) 

=GOTO(B43l) 

=RETURN0 


B.  INPUT  DIALOG  BOX 


3 

9  822 

353  Input  Sheet 

5 

500 

10 

Run  Date: 

5 

25 

10 

Item  Name: 

5 

10 

30 

Initial  Inventory  Quantities 

6 

10 

50 

On  Hand: 

5 

10 

70 

Requisition  Frequency 

6 

10 

90 

Distribution: 

6 

10 

no 

Mean- 

6 

10 

130 

Standard  Deviation; 

5 

10 

150 

Requisition  Quantity 

6 

10 

170 

Distribution: 

6 

10 

190 

Mean: 

6 

10 

210 

Standard  Deviation: 

5 

10 

230 

Stockout  Information 

5 

10 

250 

Type  Distribution: 

5 

10 

270 

Mean: 

6 

10 

290 

Standard  Deviation: 

64 


6 

10 

310 

Cost  per  Stockout: 

5 

405 

70 

"IVpe  of  Inventory  Replenishment 

5 

405 

90 

System: 

5 

405 

110 

Frequency: 

5 

215 

50 

On  Order: 

5 

419 

50 

On  Backorder: 

5 

405 

130 

Q 

5 

405 

160 

R 

5 

405 

170 

Procurement  Lead  Time 

5 

405 

190 

Type  Distribution: 

5 

405 

210 

Mean: 

5 

405 

230 

Standard  Deviation: 

5 

405 

250 

Time  Period  to  Record  Data 

6 

405 

270 

Start  Time: 

5 

406 

330 

Options: 

5 

611 

270 

End  Time: 

1 

475 

327 

88 

OK 

2 

575 

327 

88 

Cancel 

24 

675 

327 

88 

Help 

6 

135 

7  332 

Component  XYZ 

209 

585 

7  : 

190 

4/5/1993 

6 

122 

47 

60 

72 

6 

324 

47 

60 

0 

6 

567 

47 

60 

0 

21 

182 

87 

160 

rlc20:rl0c204 

8 

182 

107 

60 

0.2 

8 

182 

127 

60 

0.8 

21 

182 

167 

160 

R1C20:R10C204 

8 

182 

187 

60 

2 

8 

182 

207 

60 

6 

21 

182 

247 

160 

R1C20:R10C204 

8 

182 

267 

60 

0.8 

8 

182 

287 

60 

0.2 

8 

182 

307 

100 

100 

21 

600 

87 

160 

R1C21:R4C211 

8 

600 

107 

160 

4 

8 

600 

127 

60 

64 

8 

600 

147 

60 

18 

21 

600 

187 

160 

R1C20:R10C205 

8 

600 

207 

60 

5 

8 

600 

227 

60 

0 

8 

610 

267 

100 

1 

8 

710 

267 

100 

6 
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APPENDIX  B  AR(1)  MODEL 


A.  MAIN  MACRO 

wwnTnand  NewValues  Generates  the  bask  table 

=HIDE0 

=OPEN("working.xls‘') 

-0.9  =F0RMULA(A4;D:\THESIS\W0RKAR1  \W0RKING.XLS’!C2) 

=HIDE0 

OPENC'graphps.xIs") 

=HIDE0 

=UNHIDE(''GRAPHPS.XLS  Chart  1") 

=FOR("howM, 100,1) 

=ACTIVATE("WORKING.XLS'') 

=CALCULATE.NOW0 
=  FOR(''counter2",2,lO,l) 

=  IF(counter2=4.GOTO($B$22)) 

=  IF(coxinter2=5,GOTO($B$22)) 

=  IF(counter2=6,GOTO($B$22)) 

=  IF(counter2=8.GOTO($B$22)) 

=  IF(counter2=9,GOTO($B$22)) 

=FORMULA(GET.CELL(5,TEXTREF(”working.xls!r27c"& 

TEXT(counter2,0))),'IEXTREF("graphps.xJs!r"& 

TEXT(how+4,0)&''c"&TEXT(((counter2-l)*4)-3,0))) 

=FORMULA(GET.CELL(5,TEXTREFrworkmg.xIs!r3lc"& 

TEXT(couiiter2,0))),TEXTREF("graphps.xIs!r''& 

TEXT(how+4,0)&"c"&TEXT(((counter2-l)*4)-2,0))) 

=FORMULA(GET.CELL(6,TEXTREF("working.xls!r35c''& 

TEXT(counter2,0))),'IEXTREF("graphps.xl8!r"& 

TEXT(how+4,0)&''c''&TEXT(((co\mter2-l)*4)-l,0))) 

=FORMULA(GET.CELL(6,TEXTREF("workmg.xls!r39c"& 

TEXT(counter2,0))),TEXTREFC’graphps.xls!r"& 

TEXTChow+4,0)&"c"&TEXT(((comiter2-l)*4)-0,0))) 

=  NEXTO 

=NEXT0 

=RETURN0 
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B.  INPUT  DIALOG  CODE 


492  135 


5 

259 

3 

Parameter  Box 

5 

146 

38 

Mu  Value 

5 

127 

59 

Alpha  Value 

5 

12 

79 

Number  of  Values  required: 

5 

30 

100 

Number  of  Runs  Desired; 

8 

233 

34 

96 

0 

8 

233 

55 

96 

0 

7 

233 

76 

96 

0 

7 

233 

97 

96 

10 

1 

369 

29 

88 

OK 

2 

373 

59 

Cancel 

C.  INPUT  WORKSHEET 

Alpha  0.9000 
Mu  0.0000 


Random  16  By  By  By 


Table  As  Is 

0.728039015  0.7280 
0.293417432  0.2934 
0.702409882  0.7024 
0.360505867  0.3605 
1.146961164  1.1470 
1.137490137  1.1375 
1.109449402  1.1094 
1.471979222  1.4720 
1.347174986  1.3472 
1.486807775  1.4868 
2.163198621  2.1532 
1.576787242  1.5768 
1.228616318  1.2286 
1.527331639  1.5273 
1.5031517  1.5032 

0.958851815  0.9589 
0.899218644 


!'s  32's  256’s 

0.5107  1.0500  0.3122 
0.5315  -0.3261  0.1452 
1.1422  0.2889  0.3244 
1.2907  0.3439  0.0888 
1.4165  1.4928  -0.2324 
1.8650  -1.2448  -0.0975 
1.3780  0.2834  0.1256 
1.2310  0.6098  -0.4336 
1.1471  0.6800  -0.6756 
0.6101  1.0451  -0.0360 
0.6677  0.9943  -0.0656 
0.6526  0.9129  -0.0856 
1.2237  -0.7779  0.0710 
1.6476  0.1417  0.0877 
0.9226  -1.0406  -0.2253 
0.5633  -0.7939  0.0878 
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1.396063371  Mean 
0.451370896  1.1707 
0.768843129 
0.840207352 
0.49524308  St 

0.716012221  0.1205 
0.589230872 
1.233709939 
1.213684499  UCL 
1.581244831  1.3819 
1.713866226 
1.105883467 
0.739345047  LCL 
0.721061958  0.9696 
0.406462228 
-0.602433799 
-0.605793575 
-0.595680786 

D.  OUTPUT  GRAPHS 


Mean  Mean  Mean 
1.0500  0.2287  -0.0381 


St  Dev  St  Dev  St  Dev 
0.1063  0.2094  0.0651 


UCL  UCL  UCL 
1.2364  0.6967  0.0760 


LCL  LCL  LCL 
1.8637  -0.1383  -0.1521 
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Figure  22:  AR(1)  Model,  ♦=0.0 


69 


Figure  23:  AR(1)  Model, 


Figure  24:  AR(1)  Model,  ^-0.9 


APPENDIX  C  EAR(1)  MODEL 


A.  MAIN  MACRO 

oommand  NewVahiea  GenarateB  the  bamc  table 

=HIDEO 

=OPEN("working.xls‘') 

-0.9  =FORMULA(A4,'D:\THESIS\WORKEARl  \WORKING.XLS'!D2) 
=HIDE0 

=OPEN("  graphps.xls'') 

=HIDE0 

=UNHIDErGRAPHPS.XLS  Chart  1") 

=FOR(''howM,  100,1) 

=ACTIVATE("WORKING.XLS") 

=CALCULATE.NOW0 
=  POR("counter2'',2,lO,l) 

=  IF(counter2=4,GOTO($B$22)) 

=  IF(counter2=6,GOTO($B$22)) 

=  IF(comiter2=6,GOTO($B$22)) 

=  IF(counter2=8,GOTO($B$22)) 

*  IF(coxmter2=9,GOTO($B$22)) 

=FORMULA(GET.CELL(6,TEXTREFC'working.xls!r27c"& 
TEXT(co\mter2+ 1 ,0))),TBXTREF(''graphp8.xl8!r”& 
TEXT(how+4,0)&"c"&TEXT(((counter2-l)*4)-3,0))) 
=FORMULA(GET.CEm6,TEXTREF(’'working.xl8!r31c"& 
TEXT(coxmter2+ 1 ,0))),TBXTREFC'graphp8.xl8!r"& 
TEXT(how+4,0)&"c"&TEXT(((counter2-l)*4)-2,0))) 
=FORMULA(GET.CELL(6,TEXTREF("working.xl8!r35c"& 
TEXT(coiinter2+ 1 ,0))),TEXTREF("graphp8.xl8!r"& 
TEXT(how+4,0)&"c''&TEXT(((counter2-l)*4)-l,0))) 
=FORMULA(GET.CELL(6,TEXTREF("workiiig.xl8!r39c"& 
TEXT(counter2+l,0))),TEXTREFCgraphp8.xls!r''& 
TEXT(how+4,0)&'’c"&TEXT(((counter2-l)*4)-0,0))) 

=  NEXTO 

=NEXTO 

=RETURNO 
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B. 

INPUT  DIALOG  BOX 

492  136 

6 

259 

3 

Parameter  Box 

6 

146 

38 

Mu  Value 

5 

127 

59 

Alpha  Value 

5 

12 

79 

Number  of  Values  required: 

6 

30 

100 

Number  of  Runs  Desired: 

8 

233 

34  96 

0 

8 

233 

66  96 

0 

7 

233 

76  96 

0 

7 

233 

97  96 

10 

1 

369 

29  88 

OK 

2 

373 

69 

Cancel 

C. 

OUTPUT  GRAPHS 

ngure  25:  EAR(l)  Model,  ^-0.9 


72 


Figure  26:  EAR(l)  Model,  ♦=0.5 


Figure  27:  CAR(l)  Model,  ♦=0.0 
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APPENDIX  D  WWl  QUEUING  MODEL 


A.  MAIN  MACRO 

command  NewValues  Genearates  the  bask  table 

=DIAIX)G.BOX(ParameterBox)Calls  entry  Dialog  Box 
=HIDE0 

=OPEN("graphps.xl8  ) 

=WINDOW.MAXIMIZE0 

=HIDE0 

=UNHIDE("graphpa.xl8  Chart  I") 

=WINDOW.MAXIMIZE0 
=FORMULA("'o90''  ,C  106) 

=SET.VALUE(Arate,0.8) 

=SET.VALUE(Srate,l) 

=SET.VALUE(RHO,Arat^rate) 

=SET.VALUE(U,RANDO) 

=PORMULA(SaveCounter+ 1  ,SaveCounter) 
=OPENC'D:\the8i8\mmlpro\WORKING.XLS") 

=HIDE0 

=FORMULA(IF(U<1-RHO,0,*LN((1  -UVRHOKSrate- 

Arate)),'D:\THESIS\PART2\MMlPRO\WORKING.XLS'!A5) 
=FOR("  counter",  1  ,NumberOfV  alue8, 1) 
=SET.VALUE(Ratio,Srate(^Srate+Arate)) 

=set.valuecu,rando) 

=IF(U<Ratio) 

=  SET.VALUE(DUM,LN(U/RatioyArate) 

=ELSE0 

=  SET.VALUE(DUM,-LN((l-UKl-Ratio)ySrate) 

=END.IF0 

=FORMULA(IF(GET.CELL(5,TEXTREF("working.xl8!r"& 
TEXT(counter+4,0)&"cr))+DUM>=O.GET. 
CEm6,TEXTREF("working.xl8!r"&TEXT(counter 
+4,0)&"cl"))+DUM,0),TEXTREF("workingjd8!r"& 
TEXT(counter+ 1  +4,0)&"cl")) 

=NEXT0 

=FOR("counter"  ,1,16,1) 

=  FORMULA(TOXTREFC'working.xl8!r"&counter+4& 
+cl"),TEXTREF("working.xIs!r"&countei*+4&"c2")) 


74 


=NEXTO 

=FOR(''couiiter" ,  1 ,64 .4) 

=  FORMULA((SUM(TOXTREF(''working.xl8!r"&counter 
+4&''cl");TOXTREF("working.xls!r"&coimter 
+4+3&"cl")y4),TEXTREF(''workiiig.xl8!r"& 
((coimter+3y4)+4&"c4")) 

=NEXTO 

=FOR("counter'' , 1 ,4096,266) 

=  FORMlJLA((SUM(TEXTOEFC'working.xl8!r"&coimter+4& 
■'cl''):TEXTREF("working.xl8!r"&couiiter+4+265& 
''cl")y266),TBXTREF(''workiiig.xl8!r"&((counter 
+266V256)+4&"cl0")) 

=NEXT0 

=ACTIVATE("WORKING.XLS'') 

=SAVE.ASC'D;\the8is\mmlpro\"&Cl06&TEXT 

(SaveCouiiter,0)&''.XLSM,"",FALSE;'".FALSE) 

=HIDE0 

=  FOR("counter2",2,l0.1) 

=  IF(counter2=3,GOTO($B$67)) 

=  IF(counter2=6,GOTO($B$67))=ACTIVATE(Cl06& 
SaveCounter&”  .xls") 

=  IF(counter2=6,GOTO($B$67)) 

=  FORMULA(600+''r3lc"&TBXT(counter2,0),TEXTREF 
("graphp8.xls!r’'&TEXT(SaveCounter+4,0)&’'c”& 
TEXT(((coiinter2-l)’4)-l,0))) 

=  IF(counter2=7,GOTO($B$67))=COPY0 
=  IF(counter2=8,GOTO($B$67))=ACTIVATE(  ”GRAPHPS.X1^") 
=  IF(co\mter2=9,GOTO($B$67))=SELECT("R9C30") 
=FORMULA(600+GET.CELL(6,TEXTREF(C106&TEXT 
(SaveCounter,0)&".xl8!r23c"&TEXT(counter2, 
0))),TEXTREF("graphps.xls!r‘'&TEXT(SaveCounter 
+4,0)&"c’'&TEXT(((counter2- 1 )  •4)-3,0)))=PASTE0 
=FORMULA(GET.CELL(6,TEXTREF(Cl06&TEXTrSaveCounter, 
0)&",xl8!r27c"&TEXT(counter2,0))),TEXTREF( 
'*graphps.xls!r"&TEX'IXSaveCoimter+4,0)&’'c"& 
TEXT(((co\inter2-l)*4)-2,0)))=RETURN0 
=FORMULA(600+GET.CELL(6,TEXTREF(C106&TEXT 
(SaveCoimter,0)&".xls!r3lc”&TEXT(counter2,0))), 
TEXTREF("graphp8.xls!r’'&TEXT(SaveCounter+4, 
0)&"c"&TEXT(((couiiter2- 1)*4)- 1 ,0))) 
=FORMULA(500+GET.CELL(5,TEXTREF(C  106&TEXT 
(SaveCounter,0)&".xls!r36c.”&TEXT(coiuiter2, 
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0))),TEXTREF("graphps.xIs!r’'&:TEXT(SaveCounter 

+4,0)&''c"&TEXT(((counter2-l)*4)-0,0))) 

=  NEXTO 
=CLOSEO 

=ACTIVATE(''variance.xlm'*) 

=CHANGE.UNK(Cl06&TEXT(SaveCounter,0)&".xJs*', 

"workiiig.xls",l) 

=IF(SaveCo\mter<NtiinberOfR\ins,GOTO(Bl6),GOTO(B74)) 

=ACTIVATE("graphps.xIs'') 

=SAVE.AS(Cl06&"graph’') 

=CLOSEO 

=UNHIDE("variaiice.xIm") 

=FORMULA(0,E1) 

=IF(Cl06=''o90")=GET.CELL(6,Cl06) 

=  FORMULA('*'o60".Cl06) 

=SET.VALUE(Alpha,0.5) 

=ELSE.IF(Cl06="o60'') 

=  FOElMULAC"o00",Cl06) 

=SET.VALUB(Alpha,0) 

=ELSE.IF(Cl06=''o00'') 

=  FORMULAr'-06'',Cl06) 

*SET.VALUE(Alpha,-0.6) 

=ELSE.IF(Cl06="-06'') 

*  FORMULA("'-09",C106) 

=SET.VALUB(Alpha,-0.9) 

=ELSE.IF(Cl06=”-0.9'') 

=  GOTO(B96) 

=BND,IFO 

=GOTO(B11) 

=RETURNO 


ftiDCikm  (latheEGraphsCoQectfl  data  to  graph 

=OPENC'graphp8.xls") 

=FOR("counterM,20,l) 

OPENCnm"&TEXT(counter,0)&"jds") 

=HIDE0 

=  POR("counter2'',2,lO,l) 

=FORMULA(GET.CELL(6,TEXTREF(’’run''&TEXT(counter, 
0)&".xl8!r23c"&TEXT(comiter2,0))),TEXTREF 
("graphp8.xl8!r"&TEXT(counter+4,0)&''c"&  TEXT(((counter2' 

l)*4)-3,0))) 


76 


=FORMULA(GET.CELL(6,TEXTREF('’run"&TEXT(counter. 
0)&".xls!r27c"&TEXT(counter2,0))),TEXTREF 
("graphps.xls!r''&TEXT(counter+4,0)&"c"&  TEXT(((oounter2- 

l)'4)-2,0))) 

=F0RMULA(GET.CELL(6,TEXTREF(’'run"&TEXT 
(counter  ,0)&"  .xis!r3  lc''&TEXT(counter2,0))), 
TEXTREF("graphps.xls!r''&TEXT(counter+4,0)&''c''& 
TEXT(((counter2-l)*4)-1.0))) 

=FORMULA(GET.CELL(5,TEXTREF("nm"&TEXT(counter, 
0)&’'.xls!r35c"&TEXT(counter2,0))),TEXTREF( 
"graphps.xls!r''&TEX'IXcounter+4,0)&"c"&TEXT  (((counter2- 

l)*4)-0,0))) 

=NEXT0 

= ACTI  VATE(''run"&TEXT(counter,0)& '  .xls") 

=SAVEO 

=CLOSE0 

=NEXT0 

=RETURN0 


B.  WORKSHEET  CODE 

Alpha  0.9000 
Mu  0.0000 


Random 

Table 

-1.6213224 

-1.22069446 

-1.22776797 

-1.37944601 

-1.06692687 

-1.10964019 

-1.36727692 

-0.89182913 

-1.60931646 

-1.92031604 

-1.70964238 

-1.06421186 

-0.81630669 

-1.72862307 

-0.60097121 


16  By 
As  Is  2's 


By  By 
32's  266’s 


-1.6213  -1.4210  -0.2349  0.2909 
-1.2206  -1.3036  1.3614  0.0426 

-1.2278  -1.0833  1.6340  -0.1360 


-1.3794 

-1.1246 

-1.0669 

-1.7648 

-1.1096 

-1.3869 

-1.3673 

-1.2726 

-0.8918 

-0.4076 

-1.6093 

-0.0011 

-1.9203 

0.3013 

-1.7096 

0.2239 

-1.0642 

0.7604 

-0.8163 

0.6641 

-1.7286 

1.664 

-0.6010 

1.0992 

0.9038  0.1470 
0.0036  0.1148 
0.0013  0.2470 
-0.7909  0.1866 
-0.4414  -0.3072 
-0.1820  0.0136 
0.6717  -0.1767 
-0.0210  0.6186 
0.6693  -0.1686 
0.0683  -0.1766 
S  -0.0374  -0.0266 
0.1304  0.1686 
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-0.21416811  -0.2142  1.5037  -0.9496  0.0140 

0.05686469 

-0.06816206  Mean  Mean  Mean  Mean 

0.49801669  -1.2206  -0.2349  0.1667  0.0533 

0.10463922 
0.13928614 

0.30846461  St  Dev  St  Dev  St  Dev  St  Dev 
0.74661632  0.4627  1.1266  0.6963  0.2280 

0.75613639 
0.54237701 

0.58684934  UCL  UCL  UCL  UCL 

1.69281666  -1.0343  0.2284  0.4630  0.1471 

1.43644438 
1.02604077 

1.17246402  LCL  LCL  LCL  LCL 

1.74208696  -1.4067  -0.6982  -0.1197  -0.0405 
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APPENDIX  E  PROCESS  CONTROL 


A.  MAIN  MACRO 

=ACTIVATE("8qc.xl8") 

=FORMULA(0,'D:\THESIS\PART3\SQC.XLS'!L6) 

=FORMULA(0,'D:\THESIS\PART3\SQC.XLS’!M6) 

=FORMULA(0,'D:\THESIS\PART3\SQC.XLS'!L6) 

=FORMULA(0;D:\THESIS\PART3\SQC.XLS'!M6) 

=FOR(''a",l.GET.CELL(5,'D:\THESIS\PART3\SQC.XLS'!F5),l) 

=ACTIVATE("SQC.XLS") 

=SELECT("R8C2:R46C2","R46C2") 

=COPY0 

=PASTE.SPECIAL(3, 1  .FALSE.FALSE) 
=SELECT("R8C5:R46C5","R46C5") 

=COPY0 

=PASTE.SPECI  AL(3, 1 , FALSE, FALSE) 

=  FORMULA(a,'D:\THESIS\PART3\SQC.XLS'!H2) 

=  IF('D:\THESIS\PART3\SQC.XLS'!K6="  Accept") 
=FORMULACD:\THESIS\PART3\SQC.XLS'!L6+l,’D: 

\THESIS\PART3\SQC.XLS'!L5) 

=FORMULACD:\THESIS\PART3\SQC.XLS’!M5.’D:\THESIS 

\PART3\SQC.XLS'!M5) 

=  ELSE.IF('D:\THESIS\PART3\SQC.XLS’!K5="  Reject") 
=FORMULA('D:\THESIS\PART3\SQC.XLS’!L6, ’D:\THESIS 
\PART3\SQC.XLS'!L6) 

=FORMULA('D:\THESIS\PART3\SQC.XLS’!M5+l,’D:\THESIS 

\PART3\SQC.XLS'!M6) 

=  END.IFO 

=  IF('D;\THESIS\PART3\SQC.XLS'!K6="  Accept") 
=FORMULACD:\THESIS\PART3\SQC.XLS'!L6+l,'D:\THESIS 
\PART3\SQC.XLS'!L6) 

=FORMULACD:\THESIS\PART3\SQC.XLS’!M6, 'D:\THESIS 
\PART3\SQC.XLS'!M6) 

=  ELSE.IFCD:\THESIS\PART3\SQC.XLS'!K6="  Reject") 
=FORMULA('D;\THESIS\PART3\SQC.XLS'!L6,‘D:\THESIS 
\PART3\SQCJCLS’!L6) 

=FORMULA(’D:\THESIS\PART3\SQC.XLS'!M6+l, 'D:\THESIS 
\PART3\SQC.XLS’!M6) 

=  END.IFO 
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=ACT1VATE(''SQC.XLS) 
=SELECTC' R3C 1 7;  EU 1 C 1 7") 
=COPY0 

=SELECT("R8C2") 

=PASTEO 

=SELECT("  R3C 1 8:R4 1 C 1 8”) 
=COPY0 

=SELECT(‘'R8C5'') 

=PASTEO 

=NEXTO 

=RETURNO 


B.  WORKSHEET  LAYOUT 


alpha  = 

0 

n  = 

39 

mu  = 

10 

c  = 

1 

muprime  11.294 

Type  I 

0 

UCL  = 

12.576 

Type  II 

0 

8 

LCL  = 

7.4242 

replication 

8 

1000 

C-Process 

ooc- 

Process 

1 

12.85874 

1 

10.566325 

0 

2 

9.696326 

0 

11.010028 

0 

3 

10.86229 

0 

10.988479 

0 

4 

10.37666 

0 

9.9330193 

0 

6 

9.091763 

0 

12.887687 

1 

6 

9.968722 

0 

11.901967 

0 

7 

8.863496 

0 

10.138680 

0 

8 

9.437004 

0 

11.083608 

0 

9 

8.921389 

0 

10.998890 

0 

10 

10.76948 

0 

11.373073 

0 

11 

9.672391 

0 

12.284668 

0 

12 

10.48236 

0 

13.660179 

1 
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13 

9.790236 

0 

14 

9.771188 

0 

16 

9.169772 

0 

16 

9.302067 

0 

17 

9.733919 

0 

18 

10.90400 

0 

19 

9.556637 

0 

20 

11.06074 

0 

21 

10.00615 

0 

22 

10.62555 

0 

23 

1 

9.501144 

0 

24 

11.67813 

0 

26 

10.02455 

0 

26 

8.752327 

0 

27 

9.674731 

0 

28 

12.26949 

0 

29 

11.38882 

0 

30 

10.12019 

0 

31 

11.68462 

0 

32 

8.923697 

0 

33 

11.22563 

0 

34 

11.84096 

0 

36 

9.167223 

0 

36 

10.32481 

0 

37 

12.28669 

0 

38 

11.21761 

0 

39 

9.451619 

0 

rejection 

1 

11.245568  0 

10.306608  0 

11.854465  0 

11.671596  0 

12.037539  0 

8.8634499  0 

12.233179  0 

13.661651  1 

10.976206  0 

9.6877362  0 

10.868471  0 

9.6916219  0 

10.519138  0 

10.095880  0 

9.5511311  0 

12.792277  l 

13.311218  1 

10.668683  0 

9.6926170  0 

10.860681  0 

11.107024  0 

10.750177  0 

12.228601  0 

12.689703  1 

11.903699  0 

10.286871  0 

10.342672  0 

rejection  6 
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